EB3I n1 2025 scRNAseq
-
PRE-PROCESSING (III)
-
Filtering, cell cyle, cell doublets

1 PREAMBLE

1.1 Purpose of this session

This file describes the different steps to perform the third part of the single cell RNAseq data analysis training course for the EB3I n1 2025, covering these steps :

  • Filtering of bad cells and extreme low-expression features

  • Estimation of the cell cycle phase

  • Identification and removal of cell doublets

2 Start Rstudio

3 Warm-up

  • We set common parameters we will use throughout this session :
# setparam


## Set your project name
# WARNING : Do not just copy-paste this ! It's MY project name ! Put YOURS !!
project_name <- "ebaii_sc_teachers"


## Control if the project_name exists on the cluster
cat('PATH CHECK : ', dir.exists(paste0('/shared/projects/', project_name)))
Show output 
PATH CHECK :  TRUE
## Seed for the RNG
my_seed <- 1337L

## Minimum cells with counts to keep a feature
min_cells <- 5

4 Prepare the data structure

We will do the same as for former steps, just changing the session name.

4.1 Main directory

# maindir

## Preparing the path
TD_dir <- paste0("/shared/projects/", project_name, "/SC_TD")

## Creating the root directory
# dir.create(path = TD_dir, recursive = TRUE)

## Print the root directory on-screen
print(TD_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD"

4.2 Current session

# sessiondir

## Creating the session (Preproc.3) directory
session_dir <- paste0(TD_dir, "/03_Preproc.3")
dir.create(path = session_dir, recursive = TRUE)

## Print the session directory on-screen
print(session_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3"

4.3 Input directory

# indir

## Creating the INPUT data directory
input_dir <- paste0(session_dir, "/DATA")
dir.create(path = input_dir, recursive = TRUE)

## Print the input directory on-screen
print(input_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3/DATA"

4.4 Genelists directory

This is a directory where we will store additional information from knowledge bases about genes used to estimate the cell cycle phase of cells.

# resdir

res_dir <- paste0(TD_dir, "/Resources")
glist_dir <- paste0(res_dir, "/Genelists")

## Create the directory
dir.create(path = glist_dir, recursive = TRUE)

## Print the resources directory on-screen
print(glist_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/Resources/Genelists"

4.5 Output directory

# outdir

## Creating the OUTPUT data directory
output_dir <- paste0(session_dir, "/RESULTS")
dir.create(path = output_dir, recursive = TRUE)

## Print the output directory on-screen
print(output_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3/RESULTS"

5 Reload the Seurat Object

  • We can reload the object we saved at the former step
# dataload

## The latest Seurat object saved as RDS (name)
sobj_file <- "02_TD3A_S5_Metrics_31053.4587.RDS"

## The latest Seurat object saved as RDS (full path)
sobj_path <- paste0(TD_dir, 
                    "/02_Preproc.2/RESULTS/",
                    sobj_file)

force <- FALSE  ## To force a re-download of a Zenodo-hosted backup
local <- FALSE  ## To force a loading from a local backup

## In case of error/lost data : force a reload from a Zenodo backup repository
if(force) {
  zen_id <- "14035293"
  zen_backup_file <- paste0("https://zenodo.org/records/",
                            zen_id,
                            "/files/",
                            sobj_file)
  ## Recreate the expected path if it does not exist
  dir.create(path = dirname(sobj_path), recursive = TRUE)
  ## Download the file
  download.file(url = zen_backup_file,
                destfile = sobj_path)
}

## In case of error/lost data : force a reload from a local backup repository
if(local) {
  sobj_path <- paste0(
    "/shared/projects/2422_ebaii_n1/atelier_scrnaseq/TD/BACKUP/RDS/",
    sobj_file)
}

## Load the object
sobj <- readRDS(file = sobj_path)



Now, we finally have all the metrics and bias sources we could use, so we can actually get to the filtering step.

6 Filtering

6.1 Features

  • As already explained, the sparsity of the count matrix is high, very high.

  • To the point that there probably are some features that have so low expression level that they were only counted in very few cells. As such, they have absolutely no chance to characterize any cell type, nor harbor some variation between different cell types.

  • As such, we should discard these features


What are the actual dimensions of our expression data ?

# seudim1

## Dimensions of our Seurat object
dim(sobj)
Show output 
[1] 31053  4587


We want to remove the features that are expressed in less than 5 cells.

First, we quantify, for each feature, the number of cells with 0 count.

# nfeatdiff

## Compute the amount of 0s per feature
nFeat_zero <- sparseMatrixStats::rowCounts(
  x = SeuratObject::GetAssayData(
    object = sobj, 
    assay = "RNA", 
    layer = "counts"), 
  value = 0)

## Inversion : computing the number of cells with at least one count, per feature
nFeat_nonzero <- ncol(sobj) - nFeat_zero

## Identify those with at least 5 cells with expression
nFeat_keep <- nFeat_nonzero >= min_cells

## Quantify the selected ones
table(nFeat_keep)
Show output 
nFeat_keep
FALSE  TRUE 
18545 12508


And now we can discard features expressed in less than 5 cells :

# prefeatfiltviz

## Create a temporary object
sobj_VIZ_pre <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ_pre <- Seurat::NormalizeData(object = sobj_VIZ_pre, verbose = FALSE)
sobj_VIZ_pre <- Seurat::ScaleData(object = sobj_VIZ_pre, verbose = FALSE)
sobj_VIZ_pre <- Seurat::FindVariableFeatures(object = sobj_VIZ_pre, verbose = FALSE)
sobj_VIZ_pre <- Seurat::RunPCA(object = sobj_VIZ_pre, npcs = 21, verbose = FALSE)
sobj_VIZ_pre <- Seurat::RunUMAP(object = sobj_VIZ_pre, dims = c(1:20), verbose = FALSE)

## Cell space visualization
dp_pre <- Seurat::DimPlot(
  object = sobj_VIZ_pre, 
  reduction = 'umap') + Seurat::DarkTheme()
print(dp_pre)
Show plot 

## Remove temporary object
rm(sobj_VIZ_pre)


# featsel

## Restrict the Seurat object
sobj <- subset(x = sobj, features = rownames(sobj)[nFeat_keep])


What are the Seurat object dimensions, now ?

# seudim2

dim(sobj)
Show output 
[1] 12508  4587


We can visualize the cell space after this features filtering

# postfeatfiltviz

## Create a temporary object
sobj_VIZ_post <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ_post <- Seurat::NormalizeData(object = sobj_VIZ_post, verbose = FALSE)
sobj_VIZ_post <- Seurat::ScaleData(object = sobj_VIZ_post, verbose = FALSE)
sobj_VIZ_post <- Seurat::FindVariableFeatures(object = sobj_VIZ_post, verbose = FALSE)
sobj_VIZ_post <- Seurat::RunPCA(object = sobj_VIZ_post, npcs = 21, verbose = FALSE)
sobj_VIZ_post <- Seurat::RunUMAP(object = sobj_VIZ_post, dims = c(1:20), verbose = FALSE)

## Cell space visualization
dp_post <- Seurat::DimPlot(
  object = sobj_VIZ_post, 
  reduction = 'umap') + Seurat::DarkTheme()
print(dp_post)
Show plot 

## Remove temporary object
rm(sobj_VIZ_post)


Merging plots for ease of visualization :

# c5_pw

## Using the patchwork package to merge plots (and ggplot2 to add titles)
patchwork::wrap_plots(
  list(
    dp_pre & ggplot2::ggtitle(label = "BEFORE"),
    dp_post & ggplot2::ggtitle(label = "AFTER")), 
  nrow = 1)
Show plot


  • Question :

    # q_mincell

Do you see any difference when comparing before vs after features filtering ?


    # a_mincell

## . Not much changed
##
## . This is expected, as we removed features with almost no expression



6.2 Cells

We are now able to apply all the filtering strategies we established for each QC metric.

  • Apply the filter

    # cellfilt

## Seurat object BEFORE cell filtering
dim(sobj)

    Show output

     
    [1] 12508  4587
    ## Apply cell filtering
sobj <- subset(x = sobj, cells = colnames(sobj)[sobj$fail_qc == "pass"])

## Seurat object AFTER cell filtering
dim(sobj)

    Show output

     
    [1] 12508  4278


We can visualize the cell space since this cell filtering

# postcellfiltviz

## Create a temporary object
sobj_VIZ <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ <- Seurat::NormalizeData(object = sobj_VIZ, verbose = FALSE)
sobj_VIZ <- Seurat::ScaleData(object = sobj_VIZ, verbose = FALSE)
sobj_VIZ <- Seurat::FindVariableFeatures(object = sobj_VIZ, verbose = FALSE)
sobj_VIZ <- Seurat::RunPCA(object = sobj_VIZ, npcs = 21, verbose = FALSE)
sobj_VIZ <- Seurat::RunUMAP(object = sobj_VIZ, dims = c(1:20), verbose = FALSE)

## Cell space visualization
dp_VIZ <- Seurat::DimPlot(
  object = sobj_VIZ, 
  reduction = 'umap') + Seurat::DarkTheme()
print(dp_VIZ)
Show plot 

## Discard the temp object
rm(sobj_VIZ)


One can associate the visualization after filtering cells :

# cfilt_umaps

## Using the patchwork package to merge plots (and ggplot2 to add titles)
patchwork::wrap_plots(
  list(
    dp_post & ggplot2::ggtitle(label = "Features (count<5) filtered"),
    dp_VIZ & ggplot2::ggtitle(label = "Cells filtered (metrics)")), 
  nrow = 1)
Show plot


- Question

```{.r .question}
# q_vizfilt

Do you see any difference when comparing before vs after features filtering ?
```

<br>


```{.r .fold-hide .answer}
# q_postcellfilt

## . Not much changed as well (we did not discard many cells)
##
## . The biggest cluster structure seems more defined
```

<br>

7 Save the Seurat object

We will save our Seurat object that now contains filtered cells and features :

# saverds1

## Save our Seurat object (rich naming)
out_name <- paste0(
          output_dir, "/", paste(
            c("03", Seurat::Project(sobj), "S5", 
              "Metrics.Filtered", paste(
                dim(sobj), 
                collapse = '.'
              )
            ), collapse = "_"),
            ".RDS")

## Check
print(out_name)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3/RESULTS/03_TD3A_S5_Metrics.Filtered_12508.4278.RDS"
## Write on disk
saveRDS(object = sobj, 
        file = out_name)


8 Cell cycle scores

  • We are currently analyzing independent profiles from isolated cells, from a sample dissociation

  • As such, cells were most probably not synchronized, thus the effect of their cell cycle state on genes expression may be strong, to the point that it can bias the data (ie, mask some lower amplitude biological variation).

  • In order to assess (and maybe, remove) this bias, we have to quantify it.

  • We will perform this estimation thanks to heuristics based on knowledge : Seurat includes a method that evaluates the cell cycle phase of cells through scores for the S and G2M phases, each based on phase-specific gene signatures.

  • For this step, we will use additional gene lists from knowledge (cell cycle phase), hosted in a Zenodo respository (Id : 14037355)

8.1 Download gene lists

  • We will directly retrieve data from Zenodo to your input_dir :

    # dlzengl

## Zenodo ID
zen_id <- '14101506'

### Named files (will be used later on !)
cc_file <- "mus_musculus_Seurat_cc.genes_20191031.rds"

## Filename(s) to retrieve
toget_files <- c(cc_file)

## Folder to store retrieved files
local_folder <- glist_dir

## Use local backup ?
backup <- FALSE
if(backup) message("Using local backup !")

## Force download ?
force <- FALSE
if(force) message("Forcing (re)download !")

### Define remote folder
remote_folder <- if (backup) {
  "/shared/projects/2422_ebaii_n1/atelier_scrnaseq/TD/RESOURCES/GENELISTS/"
}  else {
  paste0("https://zenodo.org/records/", zen_id, "/files/")
}

### Reconstruct the input paths
remote_path <- paste0(remote_folder, "/", toget_files)

### Reconstruct the output paths
local_path <- paste0(local_folder, "/", toget_files)

## Retrieve files (if they don't exist), in loop
for (tg in seq_along(toget_files)) {
  ## If the file does not locally exist
  if (!file.exists(local_path[tg]) | force) {
    ## Retrieve data
    if(backup) {
      file.copy(from = remote_path[tg],
                to = local_path[tg])
    } else {
      download.file(url = remote_path[tg], 
                    destfile = local_path[tg])
    }
    ## Check if downloaded files exist locally
    if(file.exists(local_path[tg])) message("\tOK")
  } else message(paste0(toget_files[tg], " already downloaded !"))
}


8.2 Load gene lists

# cc_load

## The cell cycle gene lists file
cc_file <- paste0(glist_dir, 
                  "/", 
                  cc_file)

## Load the cell cycle reference genes lists
cc_genes <- readRDS(file = cc_file)

## Have a look on its content
str(cc_genes)
Show output 
List of 2
 $ s.genes  : chr [1:43] "Mcl5" "Pcna" "Tyms" "Fen1" ...
 $ g2m.genes: chr [1:54] "Hmgb2" "Cdk1" "Nusap1" "Ube2c" ...


  • Question

    # q_cc_check

As explained, we will use gene lists extracted from community knowledge. 
Our data contain values for genes also, so we will cross them. 

Is there something we should check ?


    # a_cc_check

## . We may check if the genes in our gene lists are effectively 
##   present in our Seurat object !
##
## . This is expected, as we removed features with almost no expression


Beyond :

  1. Write a code that performs this check

  2. Add a code to adjust the content of the gene lists accordingly (remove genes from the gene lists that are not present in our dataset)

  3. (NOTE : this is actually not needed as already checked and corrected by the cell cyle estimation method we will use)


# b_cc_answer

## Check if our data genes cover the gene lists
lapply(cc_genes, function(x) x %in% rownames(sobj))

## Remove genes not in sobj
cc_genes <- lapply(cc_genes, function(gl) { gl[gl %in% rownames(sobj)] })

## Check the modification
str(cc_genes)

## Check that all genes are available in sobj
all(unique(unlist(cc_genes)) %in% rownames(sobj))



8.3 Estimation

  • Let’s perform this estimation. But how ?

    # h_CellCycleScoring

## Reading the function help page
?Seurat::CellCycleScoring




  • Run the cell-cycle estimation

    # cc_run

## Creating a temporary object
sobj_CCS <- sobj

sobj_CCS <- Seurat::NormalizeData(object = sobj_CCS, verbose = FALSE)

## Perform the estimation
set.seed(my_seed)
sobj_CCS <- Seurat::CellCycleScoring(
  object = sobj_CCS, 
  s.features = cc_genes$s.genes, 
  g2m.features = cc_genes$g2m.genes, 
  # assay = 'RNA, ## Seurat v4
  layer = 'RNA', ## Seurat v5
  nbin = 24, 
  seed = my_seed)

## Transfering score from the temp object to the real one
sobj$CC_Seurat_S.Score <- sobj_CCS$S.Score
sobj$CC_Seurat_G2M.Score <- sobj_CCS$G2M.Score
sobj$CC_Seurat_Phase <- as.factor(sobj_CCS$Phase)

## Add "S minus G2M" score
sobj$CC_Seurat_SmG2M.Score <- sobj$CC_Seurat_S.Score - sobj$CC_Seurat_G2M.Score

## Discard temp object
rm(sobj_CCS)


  • Description of the object to see the data added

    # desc_cc

View()


8.4 Visualization

As usual, we can visualize the results as violins :

#vlncc

Seurat::VlnPlot(object = sobj,
                features = c("CC_Seurat_S.Score",
                             "CC_Seurat_G2M.Score",
                             "CC_Seurat_SmG2M.Score"))
Show plot


But it’s not that easy to interpret…

Let’s plot it in the cell space

8.4.1 Estimated cell phase :

# vizccp

## Creating a temporary object.
## This time we will keep it to speed up further plots (next chunks)
sobj_VIZ_cc <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ_cc <- Seurat::NormalizeData(object = sobj_VIZ_cc, verbose = FALSE)
sobj_VIZ_cc <- Seurat::ScaleData(object = sobj_VIZ_cc, verbose = FALSE)
sobj_VIZ_cc <- Seurat::FindVariableFeatures(object = sobj_VIZ_cc, verbose = FALSE)
sobj_VIZ_cc <- Seurat::RunPCA(object = sobj_VIZ_cc, npcs = 21, verbose = FALSE)
sobj_VIZ_cc <- Seurat::RunUMAP(object = sobj_VIZ_cc, dims = c(1:20), verbose = FALSE)

## Cell cycle phase feature plot
dp_ccp <- Seurat::DimPlot(
  object = sobj_VIZ_cc, 
  reduction = 'umap', 
  group.by = 'CC_Seurat_Phase') + Seurat::DarkTheme()
print(dp_ccp)
Show plot


8.4.2 SmG2M (“S minus G2M”) score :

# vizccs

## SmG2M feature plot
fp_cc <- Seurat::FeaturePlot(
  object = sobj_VIZ_cc, 
  reduction = 'umap',
  features = 'CC_Seurat_SmG2M.Score') + Seurat::DarkTheme()
print(fp_cc)
Show plot 

## We can discard the temp object
rm(sobj_VIZ_cc)


## {.unnumbered}

Questions

# q_cell1

Does the structure of the cells in this representation seem to have 
a link with these cell cycle phases/scores ?


# a_cell1

## It's an absolute yes !


# q_cell2

Do you think this is the result of an artifact, or
biology-related ?


# a_cell2

##
##    ¯\_(ツ)_/¯
##


9 Save the Seurat object

We will save our Seurat object that now contains the cell cycle states/scores :

# saverds2

## Save our Seurat object (rich naming)
out_name <- paste0(
          output_dir, "/", paste(
            c("04", Seurat::Project(sobj), "S5", 
              "CC", paste(
                dim(sobj), 
                collapse = '.'
              )
            ), collapse = "_"),
            ".RDS")

## Check
print(out_name)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3/RESULTS/04_TD3A_S5_CC_12508.4278.RDS"
## Write on disk
saveRDS(object = sobj, 
        file = out_name)



10 Cell doublets

We will use two different methods to detect and remove cell doublets :

  • scds (in its “hybrid” mode) : more efficient at detecting homotypic doublets

  • scDblFinder : more efficient at detecting heterotypic doublets

None of the methods accepts a SeuratObject as input, but a SingleCellExperiment object.

Hopefully :

  • Seurat has a function to perform the conversion
  • The output results can be integrated into our Seurat object with ease

10.1 Two methods

10.1.1 scds

# scds

## Fix seed
set.seed(my_seed)

## Run scds
sobj$doublet_scds.hybrid <- unname(
  scds::cxds_bcds_hybrid(
    Seurat::as.SingleCellExperiment(
      sobj, assay = "RNA"))$hybrid_score > 1)

## Contingencies
table(sobj$doublet_scds.hybrid)
Show output 

FALSE  TRUE 
 4072   206

10.1.2 scDblFinder

# scdbl

## Fix seed
set.seed(my_seed)

## Run scDblFinder (which needs another object type)
sobj$doublet_scDblFinder <- scDblFinder::scDblFinder(
  sce = Seurat::as.SingleCellExperiment(
    x = sobj, 
    assay = "RNA"),
  returnType = "table")$class == "doublet"

## Contingencies
table(sobj$doublet_scDblFinder)
Show output 

FALSE  TRUE 
 4100   178


10.2 Merge results

We merge results of the two methods

# dblmerge

## Logical union of both methods
sobj$doublet_union <- sobj$doublet_scds.hybrid | sobj$doublet_scDblFinder

## Quantify doublets
table(sobj$doublet_union)
Show output 

FALSE  TRUE 
 4035   243


We can assess tool-specific and common doublets

# dbl_types

### Singlets by default
sobj$doublet_viz <- "singlet"

### Union
sobj$doublet_viz[sobj$doublet_union] <- "both"

### scds-specific
sobj$doublet_viz[sobj$doublet_scds.hybrid & !sobj$doublet_scDblFinder] <- "scds"

### scDblFinder-specific
sobj$doublet_viz[sobj$doublet_scDblFinder & !sobj$doublet_scds.hybrid] <- "scDblFinder"

## Convert to factor
sobj$doublet_viz <- as.factor(sobj$doublet_viz)

## Contingencies
table(sobj$doublet_viz)
Show output 

       both scDblFinder        scds     singlet 
        141          37          65        4035


Beyond : Create an upset-plot for the doublet status according to the two methods used

# b_upset

## Build a list of cells tagged by one tool, the other, or both
dbl_list <-list(
  "scds" = colnames(sobj)[sobj$doublet_scds.hybrid & !sobj$doublet_scDblFinder],
  "scDblFinder" = colnames(sobj)[!sobj$doublet_scds.hybrid & sobj$doublet_scDblFinder],
  "both" = colnames(sobj)[sobj$doublet_scds.hybrid & sobj$doublet_scDblFinder])

## Draw the upset plot
UpSetR::upset(data = UpSetR::fromList(dbl_list), 
              nintersects = NA, 
              sets = rev(names(dbl_list)),
              keep.order = TRUE,
              order.by = "freq")
Show plot


Now we can remove barcodes identified as cell doublets, and visualize the cell space before and after.

10.3 Doublets filtering

10.3.1 BEFORE

# dblviz1

### Doublets viz (before removal)

## Create a temporary object
sobj_VIZ_dbl <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ_dbl <- Seurat::NormalizeData(object = sobj_VIZ_dbl, verbose = FALSE)
sobj_VIZ_dbl <- Seurat::ScaleData(object = sobj_VIZ_dbl, verbose = FALSE)
sobj_VIZ_dbl <- Seurat::FindVariableFeatures(object = sobj_VIZ_dbl, verbose = FALSE)
sobj_VIZ_dbl <- Seurat::RunPCA(object = sobj_VIZ_dbl, npcs = 21, verbose = FALSE)
sobj_VIZ_dbl <- Seurat::RunUMAP(object = sobj_VIZ_dbl, dims = c(1:20), verbose = FALSE)

## Doublets visualization
dp_dbl <- Seurat::DimPlot(
  object = sobj_VIZ_dbl, 
  reduction = 'umap',
  group.by = 'doublet_viz') + Seurat::DarkTheme()

## Plot
print(dp_dbl)
Show plot 

## Discard temp object
rm(sobj_VIZ_dbl)

10.3.2 Remove doublets

# dblrm

## Dimensions before removal
dim(sobj)
[1] 12508  4278
## Perform the removal
# sobj <- sobj[,!sobj$doublet_union]
sobj <- subset(
  x = sobj, 
  cells = colnames(sobj)[sobj$doublet_union], 
  invert = TRUE)

## Dimensions after
dim(sobj)
[1] 12508  4035

10.3.3 AFTER

# dblviz2

### Doublets viz (after removal)

## Create a temporary object
sobj_VIZ_dblf <- sobj

## Usual set of "U can't C me" functions
sobj_VIZ_dblf <- Seurat::NormalizeData(object = sobj_VIZ_dblf, verbose = FALSE)
sobj_VIZ_dblf <- Seurat::ScaleData(object = sobj_VIZ_dblf, verbose = FALSE)
sobj_VIZ_dblf <- Seurat::FindVariableFeatures(object = sobj_VIZ_dblf, verbose = FALSE)
sobj_VIZ_dblf <- Seurat::RunPCA(object = sobj_VIZ_dblf, npcs = 21, verbose = FALSE)
sobj_VIZ_dblf <- Seurat::RunUMAP(object = sobj_VIZ_dblf, dims = c(1:20), verbose = FALSE)

## Doublets-filtered visualization
dp_dblf <- Seurat::DimPlot(
  object = sobj_VIZ_dblf, 
  reduction = 'umap',
  group.by = 'doublet_viz') + Seurat::DarkTheme()

## Plot
print(dp_dblf)
Show plot 

## Discard temp object
rm(sobj_VIZ_dblf)


Merge plots for ease of use :

# dbl_umaps

## Using the patchwork package to merge plots (and ggplot2 to add titles)
patchwork::wrap_plots(
  list(
    dp_dbl & ggplot2::ggtitle(label = "Cell doublets (unfiltered)"),
    dp_dblf & ggplot2::ggtitle(label = "Cell doublets (filtered)")), 
  nrow = 1)
Show plot


Question

# q_dblcomp

What do you observe when comparing before and after the doublets filtering ?
# a_dblcomp

## . Actually not that much !
##
## . This is due to the fact that this dataset did not contain
##   a lot of heterotypic doublets
##
## . When a dataset is heavily contaminated with lots of heterotypic
##   doublets, their removal can drastically remodel the topology.


11 Save the Seurat object

We will save our Seurat object that is now filtered for doublets :

# saverds3

## Save our Seurat object (rich naming)
out_name <- paste0(
          output_dir, "/", paste(
            c("05", Seurat::Project(sobj), "S5", 
              "Doublets.Filtered", paste(
                dim(sobj), 
                collapse = '.'
              )
            ), collapse = "_"),
            ".RDS")

## Check
print(out_name)
Show output 
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/03_Preproc.3/RESULTS/05_TD3A_S5_Doublets.Filtered_12508.4035.RDS"
## Write on disk
saveRDS(object = sobj, 
        file = out_name)





12 Rsession

# rsession

utils::sessionInfo()
Show output 
R version 4.4.1 (2024-06-14)
Platform: x86_64-conda-linux-gnu
Running under: Ubuntu 22.04.5 LTS

Matrix products: default
BLAS/LAPACK: /shared/ifbstor1/software/miniconda/envs/r-4.4.1/lib/libopenblasp-r0.3.29.so;  LAPACK version 3.12.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

time zone: Europe/Paris
tzcode source: system (glibc)

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] SeuratObject_5.1.0 sp_2.2-0          

loaded via a namespace (and not attached):
  [1] RcppAnnoy_0.0.22            splines_4.4.1              
  [3] later_1.4.2                 BiocIO_1.14.0              
  [5] bitops_1.0-9                tibble_3.2.1               
  [7] polyclip_1.10-7             pROC_1.18.5                
  [9] XML_3.99-0.18               fastDummies_1.7.5          
 [11] lifecycle_1.0.4             scDblFinder_1.18.0         
 [13] edgeR_4.2.2                 globals_0.18.0             
 [15] lattice_0.22-6              MASS_7.3-65                
 [17] magrittr_2.0.3              limma_3.60.6               
 [19] plotly_4.10.4               sass_0.4.10                
 [21] rmarkdown_2.29              jquerylib_0.1.4            
 [23] yaml_2.3.10                 metapod_1.12.0             
 [25] httpuv_1.6.15               Seurat_5.3.0               
 [27] sctransform_0.4.2           spam_2.11-1                
 [29] spatstat.sparse_3.1-0       reticulate_1.42.0          
 [31] cowplot_1.1.3               pbapply_1.7-2              
 [33] RColorBrewer_1.1-3          abind_1.4-8                
 [35] zlibbioc_1.50.0             Rtsne_0.17                 
 [37] GenomicRanges_1.56.2        purrr_1.0.4                
 [39] RCurl_1.98-1.17             BiocGenerics_0.50.0        
 [41] GenomeInfoDbData_1.2.12     IRanges_2.38.1             
 [43] S4Vectors_0.42.1            ggrepel_0.9.6              
 [45] irlba_2.3.5.1               listenv_0.9.1              
 [47] spatstat.utils_3.1-4        goftest_1.2-3              
 [49] RSpectra_0.16-2             dqrng_0.4.1                
 [51] spatstat.random_3.4-1       fitdistrplus_1.2-2         
 [53] parallelly_1.45.0           DelayedMatrixStats_1.26.0  
 [55] codetools_0.2-20            DelayedArray_0.30.1        
 [57] scuttle_1.14.0              tidyselect_1.2.1           
 [59] UCSC.utils_1.0.0            farver_2.1.2               
 [61] viridis_0.6.5               ScaledMatrix_1.12.0        
 [63] matrixStats_1.5.0           stats4_4.4.1               
 [65] rmdformats_1.0.4            spatstat.explore_3.4-3     
 [67] GenomicAlignments_1.40.0    jsonlite_2.0.0             
 [69] BiocNeighbors_1.22.0        progressr_0.15.1           
 [71] scater_1.32.1               ggridges_0.5.6             
 [73] survival_3.7-0              tools_4.4.1                
 [75] ica_1.0-3                   Rcpp_1.0.14                
 [77] glue_1.8.0                  gridExtra_2.3              
 [79] SparseArray_1.4.8           xfun_0.52                  
 [81] MatrixGenerics_1.16.0       GenomeInfoDb_1.40.1        
 [83] dplyr_1.1.4                 withr_3.0.2                
 [85] fastmap_1.2.0               bluster_1.14.0             
 [87] digest_0.6.37               rsvd_1.0.5                 
 [89] R6_2.6.1                    mime_0.13                  
 [91] scattermore_1.2             tensor_1.5                 
 [93] dichromat_2.0-0.1           spatstat.data_3.1-6        
 [95] UpSetR_1.4.0                tidyr_1.3.1                
 [97] generics_0.1.4              data.table_1.17.4          
 [99] rtracklayer_1.64.0          httr_1.4.7                 
[101] htmlwidgets_1.6.4           S4Arrays_1.4.1             
[103] uwot_0.2.3                  pkgconfig_2.0.3            
[105] gtable_0.3.6                lmtest_0.9-40              
[107] SingleCellExperiment_1.26.0 XVector_0.44.0             
[109] htmltools_0.5.8.1           dotCall64_1.2              
[111] bookdown_0.39               scales_1.4.0               
[113] Biobase_2.64.0              png_0.1-8                  
[115] spatstat.univar_3.1-3       scran_1.32.0               
[117] knitr_1.50                  rstudioapi_0.17.1          
[119] rjson_0.2.23                reshape2_1.4.4             
[121] curl_6.2.3                  nlme_3.1-165               
[123] cachem_1.1.0                zoo_1.8-14                 
[125] stringr_1.5.1               KernSmooth_2.23-24         
[127] parallel_4.4.1              miniUI_0.1.2               
[129] vipor_0.4.7                 restfulr_0.0.15            
[131] ggrastr_1.0.2               pillar_1.10.2              
[133] grid_4.4.1                  vctrs_0.6.5                
[135] RANN_2.6.2                  promises_1.3.2             
[137] BiocSingular_1.20.0         beachmat_2.20.0            
[139] xtable_1.8-4                cluster_2.1.6              
[141] beeswarm_0.4.0              scds_1.20.0                
[143] evaluate_1.0.3              locfit_1.5-9.9             
[145] cli_3.6.5                   compiler_4.4.1             
[147] Rsamtools_2.20.0            rlang_1.1.6                
[149] crayon_1.5.3                future.apply_1.11.3        
[151] labeling_0.4.3              plyr_1.8.9                 
[153] ggbeeswarm_0.7.2            stringi_1.8.7              
[155] viridisLite_0.4.2           deldir_2.0-4               
[157] BiocParallel_1.38.0         Biostrings_2.72.1          
[159] lazyeval_0.2.2              spatstat.geom_3.4-1        
[161] Matrix_1.7-3                RcppHNSW_0.6.0             
[163] patchwork_1.3.0             sparseMatrixStats_1.16.0   
[165] future_1.49.0               ggplot2_3.5.2              
[167] statmod_1.5.0               shiny_1.10.0               
[169] SummarizedExperiment_1.34.0 ROCR_1.0-11                
[171] igraph_2.1.4                bslib_0.9.0                
[173] xgboost_1.7.11.1
LS0tCnRpdGxlOiAiPENFTlRFUj5FQjNJIG4xIDIwMjUgc2NSTkFzZXE8QlI+LTxCUj48Qj5QUkUtUFJPQ0VTU0lORyAoSUlJKTwvQj48QlI+LTxCUj5GaWx0ZXJpbmcsIGNlbGwgY3lsZSwgY2VsbCBkb3VibGV0czwvQ0VOVEVSPiIKZGF0ZTogIjIwMjUtMTYtMjEuMjIiCmF1dGhvcjoKICAtIG5hbWU6ICJFQjNJIG4xIHNjUk5Bc2VxIFRlYW0iCiAgLSBuYW1lOiAiQmFzdGllbiBKT0IiCiAgICBlbWFpbDogImJhc3RpZW4uam9iQGd1c3RhdmVyb3Vzc3kuZnIiCiAgLSBuYW1lOiAiTGlsaWEgWU9VTlNJIgogICAgZW1haWw6ICJsaWxpYS55b3Vuc2lAaW5zZXJtLmZyIiAKb3V0cHV0OgogIHJtZGZvcm1hdHM6OnJlYWR0aGVkb3duOgogICAgZmlnX3dpZHRoOiA4CiAgICBmaWdfaGVpZ2h0OiA2CiAgICBoaWdobGlnaHQ6IHRhbmdvICAjIyBUaGVtZSBmb3IgdGhlIGNvZGUgY2h1bmtzCiAgICBlbWJlZF9mb250czogVFJVRQogICAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlICAjIyBBZGRzIG51bWJlciB0byBoZWFkZXJzIChzZWN0aW9ucykKICAgIHRoZW1lOiBmbGF0bHkgICMjIENTUyB0aGVtZSBmb3IgdGhlIEhUTUwgcGFnZQogICAgY29sbGFwc2VkOiB0cnVlICAjIyBCeSBkZWZhdWx0LCB0aGUgVE9DIGlzIGZvbGRlZAogICAgdG9jX2RlcHRoOiAzCiAgICBzbW9vdGhfc2Nyb2xsOiB0cnVlICMjIFNtb290aCBzY3JvbGwgb2YgdGhlIEhUTUwgcGFnZQogICAgc2VsZl9jb250YWluZWQ6IHRydWUgIyMgSW5jbHVkZXMgYWxsIHBsb3RzL2ltYWdlcyB3aXRoaW4gdGhlIEhUTUwKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUgIyMgQWRkcyBhIGJ1dHRvbiB0byBkb3dubG9hZCB0aGUgUm1kCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIHRodW1ibmFpbHM6IGZhbHNlCiAgICBsaWdodGJveDogdHJ1ZQogICAgZmlnX2NhcHRpb246IGZhbHNlCiAgICBnYWxsZXJ5OiB0cnVlCiAgICB1c2VfYm9va2Rvd246IHRydWUKYWx3YXlzX2FsbG93X2h0bWw6IHRydWUgIyMgQWxsb3cgcGxhaW4gSFRNTCBjb2RlIGluIHRoZSBSbWQKZWRpdG9yX29wdGlvbnM6IAogIG1hcmtkb3duOiAKICAgIHdyYXA6IDcyCi0tLQoKPCEtLSBrbml0IHNldHVwIC0tPgoKYGBge3Iga25pdF9zZXR1cCwgZWNobyA9IEZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoCiAgZWNobyA9IFRSVUUsICAgICAgICAjIFByaW50IHRoZSBjb2RlCiAgZXZhbCA9IFRSVUUsICAgICAgICAjIFJ1biBjb21tYW5kIGxpbmVzCiAgbWVzc2FnZSA9IEZBTFNFLCAgICAjIFByaW50IG1lc3NhZ2VzCiAgcHJvbXB0ID0gRkFMU0UsICAgICAjIERvIG5vdCBkaXNwbGF5IHByb21wdAogIGNvbW1lbnQgPSBOQSwgICAgICAgIyBObyBjb21tZW50cyBvbiB0aGlzIHNlY3Rpb24KICB3YXJuaW5nID0gRkFMU0UsICAgICMgRGlzcGxheSB3YXJuaW5ncwogIHRpZHkgPSBGQUxTRSwKICBmaWcuYWxpZ249ImNlbnRlciIsIAogICMgcmVzdWx0cyA9ICdoaWRlJywKICB3aWR0aCA9IDEwMCAgICAgICAjIE51bWJlciBvZiBjaGFyYWN0ZXJzIHBlciBsaW5lCikKYGBgCgo8IS0tIENTUyB0byBjb2xvciBjaHVua3MgYW5kIG91dHB1dHMgLS0+CgpgYGB7Y3NzLCBlY2hvPUZBTFNFfQoubm90cnVuIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiBsaWdodGdyZXkgIWltcG9ydGFudDsKICBib3JkZXI6IDNweCBzb2xpZCBibGFjayAhaW1wb3J0YW50Owp9Ci5ub3RydW5vIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiBsaWdodGdyZXkgIWltcG9ydGFudDsKICBjb2xvciA6IGJsYWNrICFpbXBvcnRhbnQ7Cn0KLnF1ZXN0aW9uIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiBhcXVhbWFyaW5lICFpbXBvcnRhbnQ7CiAgY29sb3IgOiBibGFjayAhaW1wb3J0YW50OwogIGJvcmRlcjogM3B4IHNvbGlkIGxpbWVncmVlbiAhaW1wb3J0YW50Owp9Ci5xdWVzdGlvbm8gewogIGJhY2tncm91bmQtY29sb3I6IGFxdWFtYXJpbmUgIWltcG9ydGFudDsKICBjb2xvciA6IGJsYWNrICFpbXBvcnRhbnQ7Cn0KLmFuc3dlciB7CiAgYmFja2dyb3VuZC1jb2xvcjogbmF2YWpvd2hpdGUgIWltcG9ydGFudDsKICBib3JkZXI6IDNweCBzb2xpZCBicm93biAhaW1wb3J0YW50Owp9Ci5hbnN3ZXJvIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiBuYXZham93aGl0ZSAhaW1wb3J0YW50OwogIGNvbG9yIDogYmxhY2sgIWltcG9ydGFudDsKfQouYmV5b25kIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiB2aW9sZXQgIWltcG9ydGFudDsKICBib3JkZXI6IDNweCBzb2xpZCBwdXJwbGUgIWltcG9ydGFudDsKfQouYmV5b25kbyB7CiAgYmFja2dyb3VuZC1jb2xvcjogdmlvbGV0ICFpbXBvcnRhbnQ7CiAgY29sb3IgOiBibGFjayAhaW1wb3J0YW50Owp9CmBgYAoKPCEtLSBIb29rIHRvIGhhbmRsZSBjb2RlIGJsb2NrcyBvdXRwdXQgZm9sZGluZyAtLT4KCmBgYHtyIGtuaXRfaG9vaywgZWNobyA9IEZBTFNFfQpob29rcyA9IGtuaXRyOjprbml0X2hvb2tzJGdldCgpCmhvb2tfZm9sZGFibGUgPSBmdW5jdGlvbih0eXBlKSB7CiAgZm9yY2UodHlwZSkKICBmdW5jdGlvbih4LCBvcHRpb25zKSB7CiAgICByZXMgPSBob29rc1tbdHlwZV1dKHgsIG9wdGlvbnMpCiAgICAKICAgIGlmIChpc0ZBTFNFKG9wdGlvbnNbW3Bhc3RlMCgiZm9sZC4iLCB0eXBlKV1dKSkgcmV0dXJuKHJlcykKICAgIAogICAgcGFzdGUwKAogICAgICAiPGRldGFpbHM+PHN1bW1hcnk+U2hvdyAiLCB0eXBlLCAiPC9zdW1tYXJ5PlxuXG4iLAogICAgICByZXMsCiAgICAgICJcblxuPC9kZXRhaWxzPiIKICAgICkKICB9Cn0Ka25pdHI6OmtuaXRfaG9va3Mkc2V0KAogIG91dHB1dCA9IGhvb2tfZm9sZGFibGUoIm91dHB1dCIpLAogIHBsb3QgPSBob29rX2ZvbGRhYmxlKCJwbG90IikKKQpgYGAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgo8Y2VudGVyPiFbXShlYjNpX2Jhbm5lci5wbmcpPC9jZW50ZXI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBQUkVBTUJMRQoKIyMgUHVycG9zZSBvZiB0aGlzIHNlc3Npb24KClRoaXMgZmlsZSBkZXNjcmliZXMgdGhlIGRpZmZlcmVudCBzdGVwcyB0byBwZXJmb3JtIHRoZSAqKnRoaXJkIHBhcnQqKiBvZgp0aGUgc2luZ2xlIGNlbGwgUk5Bc2VxIGRhdGEgYW5hbHlzaXMgdHJhaW5pbmcgY291cnNlIGZvciB0aGUgKipFQjNJIG4xCjIwMjUqKiwgY292ZXJpbmcgdGhlc2Ugc3RlcHMgOgoKLSAgICoqRmlsdGVyaW5nKiogb2YgYmFkIGNlbGxzIGFuZCBleHRyZW1lIGxvdy1leHByZXNzaW9uIGZlYXR1cmVzCgotICAgRXN0aW1hdGlvbiBvZiB0aGUgKipjZWxsIGN5Y2xlKiogcGhhc2UKCi0gICBJZGVudGlmaWNhdGlvbiBhbmQgcmVtb3ZhbCBvZiAqKmNlbGwgZG91YmxldHMqKgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCiMgU3RhcnQgUnN0dWRpbwoKLSAgIFVzaW5nIHRoZSBbT3Blbk9uRGVtYW5kL1JzdHVkaW8gY2hlYXQKICAgIHNoZWV0XShodHRwczovL21vb2RsZS5mcmFuY2UtYmlvaW5mb3JtYXRpcXVlLmZyL3BsdWdpbmZpbGUucGhwLzE0NzUvbW9kX2ZvbGRlci9jb250ZW50LzAvT29EX1JfUnN0dWRpby5odG1sKXt0YXJnZXQ9Il9ibGFuayJ9LAogICAgY29ubmVjdCB0byB0aGUgW09wZW5PbkRlbWFuZAogICAgcG9ydGFsXShodHRwczovL29uZGVtYW5kLmNsdXN0ZXIuZnJhbmNlLWJpb2luZm9ybWF0aXF1ZS5mcil7dGFyZ2V0PSJfYmxhbmsifSBhbmQKICAgICoqY3JlYXRlIGEgUnN0dWRpbyBzZXNzaW9uKiogd2l0aCB0aGUgcmlnaHQgcmVzb3VyY2UgcmVxdWlyZW1lbnRzLgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCiMgV2FybS11cAoKLSAgIFdlIHNldCAqKmNvbW1vbiBwYXJhbWV0ZXJzKiogd2Ugd2lsbCB1c2UgdGhyb3VnaG91dCB0aGlzIHNlc3Npb24gOgoKYGBge3Igc2V0cGFyYW19CiMgc2V0cGFyYW0KCgojIyBTZXQgeW91ciBwcm9qZWN0IG5hbWUKIyBXQVJOSU5HIDogRG8gbm90IGp1c3QgY29weS1wYXN0ZSB0aGlzICEgSXQncyBNWSBwcm9qZWN0IG5hbWUgISBQdXQgWU9VUlMgISEKcHJvamVjdF9uYW1lIDwtICJlYmFpaV9zY190ZWFjaGVycyIKCgojIyBDb250cm9sIGlmIHRoZSBwcm9qZWN0X25hbWUgZXhpc3RzIG9uIHRoZSBjbHVzdGVyCmNhdCgnUEFUSCBDSEVDSyA6ICcsIGRpci5leGlzdHMocGFzdGUwKCcvc2hhcmVkL3Byb2plY3RzLycsIHByb2plY3RfbmFtZSkpKQoKIyMgU2VlZCBmb3IgdGhlIFJORwpteV9zZWVkIDwtIDEzMzdMCgojIyBNaW5pbXVtIGNlbGxzIHdpdGggY291bnRzIHRvIGtlZXAgYSBmZWF0dXJlCm1pbl9jZWxscyA8LSA1CgpgYGAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojIFByZXBhcmUgdGhlIGRhdGEgc3RydWN0dXJlCgpXZSB3aWxsIGRvIHRoZSBzYW1lIGFzIGZvciBmb3JtZXIgc3RlcHMsIGp1c3QgY2hhbmdpbmcgdGhlIHNlc3Npb24gbmFtZS4KCiMjIE1haW4gZGlyZWN0b3J5CgpgYGB7ciBtYWluZGlyLCBmb2xkLm91dHB1dCA9IEZBTFNFfQojIG1haW5kaXIKCiMjIFByZXBhcmluZyB0aGUgcGF0aApURF9kaXIgPC0gcGFzdGUwKCIvc2hhcmVkL3Byb2plY3RzLyIsIHByb2plY3RfbmFtZSwgIi9TQ19URCIpCgojIyBDcmVhdGluZyB0aGUgcm9vdCBkaXJlY3RvcnkKIyBkaXIuY3JlYXRlKHBhdGggPSBURF9kaXIsIHJlY3Vyc2l2ZSA9IFRSVUUpCgojIyBQcmludCB0aGUgcm9vdCBkaXJlY3Rvcnkgb24tc2NyZWVuCnByaW50KFREX2RpcikKCmBgYAoKIyMgQ3VycmVudCBzZXNzaW9uCgpgYGB7ciBzZXNzaW9uZGlyLCBmb2xkLm91dHB1dCA9IEZBTFNFfQojIHNlc3Npb25kaXIKCiMjIENyZWF0aW5nIHRoZSBzZXNzaW9uIChQcmVwcm9jLjMpIGRpcmVjdG9yeQpzZXNzaW9uX2RpciA8LSBwYXN0ZTAoVERfZGlyLCAiLzAzX1ByZXByb2MuMyIpCmRpci5jcmVhdGUocGF0aCA9IHNlc3Npb25fZGlyLCByZWN1cnNpdmUgPSBUUlVFKQoKIyMgUHJpbnQgdGhlIHNlc3Npb24gZGlyZWN0b3J5IG9uLXNjcmVlbgpwcmludChzZXNzaW9uX2RpcikKCmBgYAoKIyMgSW5wdXQgZGlyZWN0b3J5CgpgYGB7ciBpbmRpciwgZm9sZC5vdXRwdXQgPSBGQUxTRX0KIyBpbmRpcgoKIyMgQ3JlYXRpbmcgdGhlIElOUFVUIGRhdGEgZGlyZWN0b3J5CmlucHV0X2RpciA8LSBwYXN0ZTAoc2Vzc2lvbl9kaXIsICIvREFUQSIpCmRpci5jcmVhdGUocGF0aCA9IGlucHV0X2RpciwgcmVjdXJzaXZlID0gVFJVRSkKCiMjIFByaW50IHRoZSBpbnB1dCBkaXJlY3Rvcnkgb24tc2NyZWVuCnByaW50KGlucHV0X2RpcikKCmBgYAoKIyMgR2VuZWxpc3RzIGRpcmVjdG9yeQoKVGhpcyBpcyBhIGRpcmVjdG9yeSB3aGVyZSB3ZSB3aWxsIHN0b3JlIGFkZGl0aW9uYWwgaW5mb3JtYXRpb24gZnJvbSAKa25vd2xlZGdlIGJhc2VzIGFib3V0IGdlbmVzIHVzZWQgdG8gZXN0aW1hdGUgdGhlIGNlbGwgY3ljbGUgcGhhc2Ugb2YgY2VsbHMuCgpgYGB7ciByZXNkaXIsIGZvbGQub3V0cHV0ID0gRkFMU0V9CiMgcmVzZGlyCgpyZXNfZGlyIDwtIHBhc3RlMChURF9kaXIsICIvUmVzb3VyY2VzIikKZ2xpc3RfZGlyIDwtIHBhc3RlMChyZXNfZGlyLCAiL0dlbmVsaXN0cyIpCgojIyBDcmVhdGUgdGhlIGRpcmVjdG9yeQpkaXIuY3JlYXRlKHBhdGggPSBnbGlzdF9kaXIsIHJlY3Vyc2l2ZSA9IFRSVUUpCgojIyBQcmludCB0aGUgcmVzb3VyY2VzIGRpcmVjdG9yeSBvbi1zY3JlZW4KcHJpbnQoZ2xpc3RfZGlyKQoKYGBgCgojIyBPdXRwdXQgZGlyZWN0b3J5CgpgYGB7ciBvdXRkaXIsIGZvbGQub3V0cHV0ID0gRkFMU0V9CiMgb3V0ZGlyCgojIyBDcmVhdGluZyB0aGUgT1VUUFVUIGRhdGEgZGlyZWN0b3J5Cm91dHB1dF9kaXIgPC0gcGFzdGUwKHNlc3Npb25fZGlyLCAiL1JFU1VMVFMiKQpkaXIuY3JlYXRlKHBhdGggPSBvdXRwdXRfZGlyLCByZWN1cnNpdmUgPSBUUlVFKQoKIyMgUHJpbnQgdGhlIG91dHB1dCBkaXJlY3Rvcnkgb24tc2NyZWVuCnByaW50KG91dHB1dF9kaXIpCgpgYGAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojIFJlbG9hZCB0aGUgU2V1cmF0IE9iamVjdAoKLSAgIFdlIGNhbiByZWxvYWQgdGhlIG9iamVjdCB3ZSBzYXZlZCBhdCB0aGUgZm9ybWVyIHN0ZXAKCmBgYHtyIGRhdGFsb2FkfQojIGRhdGFsb2FkCgojIyBUaGUgbGF0ZXN0IFNldXJhdCBvYmplY3Qgc2F2ZWQgYXMgUkRTIChuYW1lKQpzb2JqX2ZpbGUgPC0gIjAyX1REM0FfUzVfTWV0cmljc18zMTA1My40NTg3LlJEUyIKCiMjIFRoZSBsYXRlc3QgU2V1cmF0IG9iamVjdCBzYXZlZCBhcyBSRFMgKGZ1bGwgcGF0aCkKc29ial9wYXRoIDwtIHBhc3RlMChURF9kaXIsIAogICAgICAgICAgICAgICAgICAgICIvMDJfUHJlcHJvYy4yL1JFU1VMVFMvIiwKICAgICAgICAgICAgICAgICAgICBzb2JqX2ZpbGUpCgpmb3JjZSA8LSBGQUxTRSAgIyMgVG8gZm9yY2UgYSByZS1kb3dubG9hZCBvZiBhIFplbm9kby1ob3N0ZWQgYmFja3VwCmxvY2FsIDwtIEZBTFNFICAjIyBUbyBmb3JjZSBhIGxvYWRpbmcgZnJvbSBhIGxvY2FsIGJhY2t1cAoKIyMgSW4gY2FzZSBvZiBlcnJvci9sb3N0IGRhdGEgOiBmb3JjZSBhIHJlbG9hZCBmcm9tIGEgWmVub2RvIGJhY2t1cCByZXBvc2l0b3J5CmlmKGZvcmNlKSB7CiAgemVuX2lkIDwtICIxNDAzNTI5MyIKICB6ZW5fYmFja3VwX2ZpbGUgPC0gcGFzdGUwKCJodHRwczovL3plbm9kby5vcmcvcmVjb3Jkcy8iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgemVuX2lkLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIi9maWxlcy8iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgc29ial9maWxlKQogICMjIFJlY3JlYXRlIHRoZSBleHBlY3RlZCBwYXRoIGlmIGl0IGRvZXMgbm90IGV4aXN0CiAgZGlyLmNyZWF0ZShwYXRoID0gZGlybmFtZShzb2JqX3BhdGgpLCByZWN1cnNpdmUgPSBUUlVFKQogICMjIERvd25sb2FkIHRoZSBmaWxlCiAgZG93bmxvYWQuZmlsZSh1cmwgPSB6ZW5fYmFja3VwX2ZpbGUsCiAgICAgICAgICAgICAgICBkZXN0ZmlsZSA9IHNvYmpfcGF0aCkKfQoKIyMgSW4gY2FzZSBvZiBlcnJvci9sb3N0IGRhdGEgOiBmb3JjZSBhIHJlbG9hZCBmcm9tIGEgbG9jYWwgYmFja3VwIHJlcG9zaXRvcnkKaWYobG9jYWwpIHsKICBzb2JqX3BhdGggPC0gcGFzdGUwKAogICAgIi9zaGFyZWQvcHJvamVjdHMvMjQyMl9lYmFpaV9uMS9hdGVsaWVyX3Njcm5hc2VxL1REL0JBQ0tVUC9SRFMvIiwKICAgIHNvYmpfZmlsZSkKfQoKIyMgTG9hZCB0aGUgb2JqZWN0CnNvYmogPC0gcmVhZFJEUyhmaWxlID0gc29ial9wYXRoKQoKYGBgCgo8YnI+PGJyPgoKTm93LCB3ZSBmaW5hbGx5IGhhdmUgYWxsIHRoZSBtZXRyaWNzIGFuZCBiaWFzIHNvdXJjZXMgd2UgY291bGQKdXNlLCBzbyB3ZSBjYW4gYWN0dWFsbHkgZ2V0IHRvIHRoZSBmaWx0ZXJpbmcgc3RlcC4KCgojIEZpbHRlcmluZwoKIyMgRmVhdHVyZXMKCi0gICBBcyBhbHJlYWR5IGV4cGxhaW5lZCwgdGhlICoqc3BhcnNpdHkqKiBvZiB0aGUgY291bnQgbWF0cml4IGlzIGhpZ2gsIF92ZXJ5XwpoaWdoLiAKCi0gICBUbyB0aGUgcG9pbnQgdGhhdCB0aGVyZSBwcm9iYWJseSBhcmUgc29tZSBmZWF0dXJlcyB0aGF0IGhhdmUgKipzbwpsb3cgZXhwcmVzc2lvbioqIGxldmVsIHRoYXQgdGhleSB3ZXJlIG9ubHkgY291bnRlZCBpbiAqKnZlcnkgZmV3IGNlbGxzKiouIEFzCnN1Y2gsIHRoZXkgaGF2ZSBhYnNvbHV0ZWx5ICoqbm8gY2hhbmNlKiogdG8gY2hhcmFjdGVyaXplIGFueSBjZWxsIHR5cGUsIG5vcgpoYXJib3Igc29tZSB2YXJpYXRpb24gYmV0d2VlbiBkaWZmZXJlbnQgY2VsbCB0eXBlcy4KCi0gICBBcyBzdWNoLCB3ZSBzaG91bGQgKipkaXNjYXJkKiogdGhlc2UgZmVhdHVyZXMKCjxicj4KCldoYXQgYXJlIHRoZSBhY3R1YWwgZGltZW5zaW9ucyBvZiBvdXIgZXhwcmVzc2lvbiBkYXRhID8KCmBgYHtyIHNldWRpbTF9CiMgc2V1ZGltMQoKIyMgRGltZW5zaW9ucyBvZiBvdXIgU2V1cmF0IG9iamVjdApkaW0oc29iaikKCmBgYAoKPGJyPgoKV2Ugd2FudCB0byByZW1vdmUgdGhlIGZlYXR1cmVzIHRoYXQgYXJlIGV4cHJlc3NlZCBpbiBsZXNzIHRoYW4gKipgciBtaW5fY2VsbHNgKioKY2VsbHMuCgpGaXJzdCwgd2UgcXVhbnRpZnksIGZvciBlYWNoIGZlYXR1cmUsIHRoZSBudW1iZXIgb2YgY2VsbHMgd2l0aCBgMGAgY291bnQuCgpgYGB7ciBuZmVhdGRpZmZ9CiMgbmZlYXRkaWZmCgojIyBDb21wdXRlIHRoZSBhbW91bnQgb2YgMHMgcGVyIGZlYXR1cmUKbkZlYXRfemVybyA8LSBzcGFyc2VNYXRyaXhTdGF0czo6cm93Q291bnRzKAogIHggPSBTZXVyYXRPYmplY3Q6OkdldEFzc2F5RGF0YSgKICAgIG9iamVjdCA9IHNvYmosIAogICAgYXNzYXkgPSAiUk5BIiwgCiAgICBsYXllciA9ICJjb3VudHMiKSwgCiAgdmFsdWUgPSAwKQoKIyMgSW52ZXJzaW9uIDogY29tcHV0aW5nIHRoZSBudW1iZXIgb2YgY2VsbHMgd2l0aCBhdCBsZWFzdCBvbmUgY291bnQsIHBlciBmZWF0dXJlCm5GZWF0X25vbnplcm8gPC0gbmNvbChzb2JqKSAtIG5GZWF0X3plcm8KCiMjIElkZW50aWZ5IHRob3NlIHdpdGggYXQgbGVhc3QgNSBjZWxscyB3aXRoIGV4cHJlc3Npb24KbkZlYXRfa2VlcCA8LSBuRmVhdF9ub256ZXJvID49IG1pbl9jZWxscwoKIyMgUXVhbnRpZnkgdGhlIHNlbGVjdGVkIG9uZXMKdGFibGUobkZlYXRfa2VlcCkKCmBgYAoKPGJyPgoKQW5kIG5vdyB3ZSBjYW4gKipkaXNjYXJkKiogZmVhdHVyZXMgZXhwcmVzc2VkIGluIGxlc3MgdGhhbiBgciBtaW5fY2VsbHNgIGNlbGxzIDoKCmBgYHtyIHByZWZlYXRmaWx0dml6LCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgcHJlZmVhdGZpbHR2aXoKCiMjIENyZWF0ZSBhIHRlbXBvcmFyeSBvYmplY3QKc29ial9WSVpfcHJlIDwtIHNvYmoKCiMjIFVzdWFsIHNldCBvZiAiVSBjYW4ndCBDIG1lIiBmdW5jdGlvbnMKc29ial9WSVpfcHJlIDwtIFNldXJhdDo6Tm9ybWFsaXplRGF0YShvYmplY3QgPSBzb2JqX1ZJWl9wcmUsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfcHJlIDwtIFNldXJhdDo6U2NhbGVEYXRhKG9iamVjdCA9IHNvYmpfVklaX3ByZSwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9wcmUgPC0gU2V1cmF0OjpGaW5kVmFyaWFibGVGZWF0dXJlcyhvYmplY3QgPSBzb2JqX1ZJWl9wcmUsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfcHJlIDwtIFNldXJhdDo6UnVuUENBKG9iamVjdCA9IHNvYmpfVklaX3ByZSwgbnBjcyA9IDIxLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX3ByZSA8LSBTZXVyYXQ6OlJ1blVNQVAob2JqZWN0ID0gc29ial9WSVpfcHJlLCBkaW1zID0gYygxOjIwKSwgdmVyYm9zZSA9IEZBTFNFKQoKIyMgQ2VsbCBzcGFjZSB2aXN1YWxpemF0aW9uCmRwX3ByZSA8LSBTZXVyYXQ6OkRpbVBsb3QoCiAgb2JqZWN0ID0gc29ial9WSVpfcHJlLCAKICByZWR1Y3Rpb24gPSAndW1hcCcpICsgU2V1cmF0OjpEYXJrVGhlbWUoKQpwcmludChkcF9wcmUpCgojIyBSZW1vdmUgdGVtcG9yYXJ5IG9iamVjdApybShzb2JqX1ZJWl9wcmUpCgpgYGAKCjxicj4KCmBgYHtyIGZlYXRzZWx9CiMgZmVhdHNlbAoKIyMgUmVzdHJpY3QgdGhlIFNldXJhdCBvYmplY3QKc29iaiA8LSBzdWJzZXQoeCA9IHNvYmosIGZlYXR1cmVzID0gcm93bmFtZXMoc29iailbbkZlYXRfa2VlcF0pCgpgYGAKCjxicj4KCldoYXQgYXJlIHRoZSBTZXVyYXQgb2JqZWN0IGRpbWVuc2lvbnMsIG5vdyA/CgpgYGB7ciBzZXVkaW0yLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgc2V1ZGltMgoKZGltKHNvYmopCgpgYGAKCjxicj4KCldlIGNhbiAqKnZpc3VhbGl6ZSB0aGUgY2VsbCBzcGFjZSoqIGFmdGVyIHRoaXMgZmVhdHVyZXMgZmlsdGVyaW5nCgpgYGB7ciBwb3N0ZmVhdGZpbHR2aXp9CiMgcG9zdGZlYXRmaWx0dml6CgojIyBDcmVhdGUgYSB0ZW1wb3Jhcnkgb2JqZWN0CnNvYmpfVklaX3Bvc3QgPC0gc29iagoKIyMgVXN1YWwgc2V0IG9mICJVIGNhbid0IEMgbWUiIGZ1bmN0aW9ucwpzb2JqX1ZJWl9wb3N0IDwtIFNldXJhdDo6Tm9ybWFsaXplRGF0YShvYmplY3QgPSBzb2JqX1ZJWl9wb3N0LCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX3Bvc3QgPC0gU2V1cmF0OjpTY2FsZURhdGEob2JqZWN0ID0gc29ial9WSVpfcG9zdCwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9wb3N0IDwtIFNldXJhdDo6RmluZFZhcmlhYmxlRmVhdHVyZXMob2JqZWN0ID0gc29ial9WSVpfcG9zdCwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9wb3N0IDwtIFNldXJhdDo6UnVuUENBKG9iamVjdCA9IHNvYmpfVklaX3Bvc3QsIG5wY3MgPSAyMSwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9wb3N0IDwtIFNldXJhdDo6UnVuVU1BUChvYmplY3QgPSBzb2JqX1ZJWl9wb3N0LCBkaW1zID0gYygxOjIwKSwgdmVyYm9zZSA9IEZBTFNFKQoKIyMgQ2VsbCBzcGFjZSB2aXN1YWxpemF0aW9uCmRwX3Bvc3QgPC0gU2V1cmF0OjpEaW1QbG90KAogIG9iamVjdCA9IHNvYmpfVklaX3Bvc3QsIAogIHJlZHVjdGlvbiA9ICd1bWFwJykgKyBTZXVyYXQ6OkRhcmtUaGVtZSgpCnByaW50KGRwX3Bvc3QpCgojIyBSZW1vdmUgdGVtcG9yYXJ5IG9iamVjdApybShzb2JqX1ZJWl9wb3N0KQoKYGBgCgo8YnI+CgpNZXJnaW5nIHBsb3RzIGZvciBlYXNlIG9mIHZpc3VhbGl6YXRpb24gOgoKYGBge3IgYzVfcHcsIGZpZy53aWR0aCA9IDE1LCBmaWcuaGVpZ2h0ID0gNiwgY2xhc3Muc291cmNlPWMoImZvbGQtaGlkZSIsICJub3RydW4iKSwgY2xhc3Mub3V0cHV0PSJub3RydW5vIn0KIyBjNV9wdwoKIyMgVXNpbmcgdGhlIHBhdGNod29yayBwYWNrYWdlIHRvIG1lcmdlIHBsb3RzIChhbmQgZ2dwbG90MiB0byBhZGQgdGl0bGVzKQpwYXRjaHdvcms6OndyYXBfcGxvdHMoCiAgbGlzdCgKICAgIGRwX3ByZSAmIGdncGxvdDI6OmdndGl0bGUobGFiZWwgPSAiQkVGT1JFIiksCiAgICBkcF9wb3N0ICYgZ2dwbG90Mjo6Z2d0aXRsZShsYWJlbCA9ICJBRlRFUiIpKSwgCiAgbnJvdyA9IDEpCgpgYGAKCjxicj4KCi0gICAqKlF1ZXN0aW9uKiogOiAKICAgIAogICAgYGBge3IgcV9taW5jZWxsLCBldmFsID0gRkFMU0UsIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiAgICAjIHFfbWluY2VsbAogICAgCiAgICBEbyB5b3Ugc2VlIGFueSBkaWZmZXJlbmNlIHdoZW4gY29tcGFyaW5nIGJlZm9yZSB2cyBhZnRlciBmZWF0dXJlcyBmaWx0ZXJpbmcgPwogICAgCiAgICBgYGAKICAgIAogICAgPGJyPgogICAgCiAgICBgYGB7ciBhX21pbmNlbGwsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX21pbmNlbGwKICAgIAogICAgIyMgLiBOb3QgbXVjaCBjaGFuZ2VkCiAgICAjIwogICAgIyMgLiBUaGlzIGlzIGV4cGVjdGVkLCBhcyB3ZSByZW1vdmVkIGZlYXR1cmVzIHdpdGggYWxtb3N0IG5vIGV4cHJlc3Npb24KICAgIGBgYAogICAgCjxicj48YnI+CgojIyBDZWxscwoKV2UgYXJlIG5vdyBhYmxlIHRvICoqYXBwbHkgYWxsIHRoZSBmaWx0ZXJpbmcgc3RyYXRlZ2llcyoqIHdlIGVzdGFibGlzaGVkIGZvcgplYWNoIFFDIG1ldHJpYy4KCgotICAgQXBwbHkgdGhlIGZpbHRlcgoKICAgIGBgYHtyIGNlbGxmaWx0fQogICAgIyBjZWxsZmlsdAogICAgCiAgICAjIyBTZXVyYXQgb2JqZWN0IEJFRk9SRSBjZWxsIGZpbHRlcmluZwogICAgZGltKHNvYmopCiAgICAKICAgICMjIEFwcGx5IGNlbGwgZmlsdGVyaW5nCiAgICBzb2JqIDwtIHN1YnNldCh4ID0gc29iaiwgY2VsbHMgPSBjb2xuYW1lcyhzb2JqKVtzb2JqJGZhaWxfcWMgPT0gInBhc3MiXSkKICAgIAogICAgIyMgU2V1cmF0IG9iamVjdCBBRlRFUiBjZWxsIGZpbHRlcmluZwogICAgZGltKHNvYmopCiAgICAKICAgIGBgYAoKPGJyPgoKV2UgY2FuICoqdmlzdWFsaXplIHRoZSBjZWxsIHNwYWNlKiogc2luY2UgdGhpcyBjZWxsIGZpbHRlcmluZwoKYGBge3IgcG9zdGNlbGxmaWx0dml6fQojIHBvc3RjZWxsZmlsdHZpegoKIyMgQ3JlYXRlIGEgdGVtcG9yYXJ5IG9iamVjdApzb2JqX1ZJWiA8LSBzb2JqCgojIyBVc3VhbCBzZXQgb2YgIlUgY2FuJ3QgQyBtZSIgZnVuY3Rpb25zCnNvYmpfVklaIDwtIFNldXJhdDo6Tm9ybWFsaXplRGF0YShvYmplY3QgPSBzb2JqX1ZJWiwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWiA8LSBTZXVyYXQ6OlNjYWxlRGF0YShvYmplY3QgPSBzb2JqX1ZJWiwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWiA8LSBTZXVyYXQ6OkZpbmRWYXJpYWJsZUZlYXR1cmVzKG9iamVjdCA9IHNvYmpfVklaLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaIDwtIFNldXJhdDo6UnVuUENBKG9iamVjdCA9IHNvYmpfVklaLCBucGNzID0gMjEsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVogPC0gU2V1cmF0OjpSdW5VTUFQKG9iamVjdCA9IHNvYmpfVklaLCBkaW1zID0gYygxOjIwKSwgdmVyYm9zZSA9IEZBTFNFKQoKIyMgQ2VsbCBzcGFjZSB2aXN1YWxpemF0aW9uCmRwX1ZJWiA8LSBTZXVyYXQ6OkRpbVBsb3QoCiAgb2JqZWN0ID0gc29ial9WSVosIAogIHJlZHVjdGlvbiA9ICd1bWFwJykgKyBTZXVyYXQ6OkRhcmtUaGVtZSgpCnByaW50KGRwX1ZJWikKCiMjIERpc2NhcmQgdGhlIHRlbXAgb2JqZWN0CnJtKHNvYmpfVklaKQoKYGBgCgo8YnI+CgpPbmUgY2FuIGFzc29jaWF0ZSB0aGUgdmlzdWFsaXphdGlvbiBhZnRlciBmaWx0ZXJpbmcgY2VsbHMgOgoKYGBge3IgY2ZpbHRfdW1hcHMsIGZpZy53aWR0aCA9IDEyLCBmaWcuaGVpZ2h0ID0gNiwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojIGNmaWx0X3VtYXBzCgojIyBVc2luZyB0aGUgcGF0Y2h3b3JrIHBhY2thZ2UgdG8gbWVyZ2UgcGxvdHMgKGFuZCBnZ3Bsb3QyIHRvIGFkZCB0aXRsZXMpCnBhdGNod29yazo6d3JhcF9wbG90cygKICBsaXN0KAogICAgZHBfcG9zdCAmIGdncGxvdDI6OmdndGl0bGUobGFiZWwgPSAiRmVhdHVyZXMgKGNvdW50PDUpIGZpbHRlcmVkIiksCiAgICBkcF9WSVogJiBnZ3Bsb3QyOjpnZ3RpdGxlKGxhYmVsID0gIkNlbGxzIGZpbHRlcmVkIChtZXRyaWNzKSIpKSwgCiAgbnJvdyA9IDEpCgpgYGAKCjxicj4KLSAgICoqUXVlc3Rpb24qKgogICAgCiAgICBgYGB7ciBxX3ZpemZpbHQsIGV2YWwgPSBGQUxTRSwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV92aXpmaWx0CiAgICAKICAgIERvIHlvdSBzZWUgYW55IGRpZmZlcmVuY2Ugd2hlbiBjb21wYXJpbmcgYmVmb3JlIHZzIGFmdGVyIGZlYXR1cmVzIGZpbHRlcmluZyA/CiAgICAKICAgIGBgYAogICAgCiAgICA8YnI+CiAgICAKICAgIGBgYHtyIHFfcG9zdGNlbGxmaWx0LCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYW5zd2VyIiksIGV2YWwgPSBGQUxTRX0KICAgICMgcV9wb3N0Y2VsbGZpbHQKICAgIAogICAgIyMgLiBOb3QgbXVjaCBjaGFuZ2VkIGFzIHdlbGwgKHdlIGRpZCBub3QgZGlzY2FyZCBtYW55IGNlbGxzKQogICAgIyMKICAgICMjIC4gVGhlIGJpZ2dlc3QgY2x1c3RlciBzdHJ1Y3R1cmUgc2VlbXMgbW9yZSBkZWZpbmVkCiAgICAKICAgIGBgYAogICAgCiAgICA8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBTYXZlIHRoZSBTZXVyYXQgb2JqZWN0CgpXZSB3aWxsICoqc2F2ZSBvdXIgU2V1cmF0IG9iamVjdCoqIHRoYXQgbm93IGNvbnRhaW5zICoqZmlsdGVyZWQgY2VsbHMgYW5kCmZlYXR1cmVzKiogOgoKYGBge3Igc2F2ZXJkczEsIGZvbGQub3V0cHV0ID0gRkFMU0V9CiMgc2F2ZXJkczEKCiMjIFNhdmUgb3VyIFNldXJhdCBvYmplY3QgKHJpY2ggbmFtaW5nKQpvdXRfbmFtZSA8LSBwYXN0ZTAoCiAgICAgICAgICBvdXRwdXRfZGlyLCAiLyIsIHBhc3RlKAogICAgICAgICAgICBjKCIwMyIsIFNldXJhdDo6UHJvamVjdChzb2JqKSwgIlM1IiwgCiAgICAgICAgICAgICAgIk1ldHJpY3MuRmlsdGVyZWQiLCBwYXN0ZSgKICAgICAgICAgICAgICAgIGRpbShzb2JqKSwgCiAgICAgICAgICAgICAgICBjb2xsYXBzZSA9ICcuJwogICAgICAgICAgICAgICkKICAgICAgICAgICAgKSwgY29sbGFwc2UgPSAiXyIpLAogICAgICAgICAgICAiLlJEUyIpCgojIyBDaGVjawpwcmludChvdXRfbmFtZSkKCiMjIFdyaXRlIG9uIGRpc2sKc2F2ZVJEUyhvYmplY3QgPSBzb2JqLCAKICAgICAgICBmaWxlID0gb3V0X25hbWUpCgpgYGAKPGJyPgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCiMgQ2VsbCBjeWNsZSBzY29yZXMKCi0gICBXZSBhcmUgY3VycmVudGx5IGFuYWx5emluZyAqKmluZGVwZW5kZW50KiogcHJvZmlsZXMgZnJvbSAgKippc29sYXRlZCoqIGNlbGxzLCBmcm9tIGEgc2FtcGxlIGRpc3NvY2lhdGlvbgoKLSAgIEFzIHN1Y2gsIGNlbGxzIHdlcmUgbW9zdCBwcm9iYWJseSAqKm5vdCBzeW5jaHJvbml6ZWQqKiwgdGh1cyB0aGUgZWZmZWN0IG9mIHRoZWlyICoqY2VsbCBjeWNsZSBzdGF0ZSoqIG9uIGdlbmVzIGV4cHJlc3Npb24gbWF5IGJlIHN0cm9uZywgdG8gdGhlIHBvaW50IHRoYXQgaXQgY2FuICoqYmlhcyB0aGUgZGF0YSoqIChpZSwgbWFzayBzb21lIGxvd2VyIGFtcGxpdHVkZSBiaW9sb2dpY2FsIHZhcmlhdGlvbikuCgotICAgSW4gb3JkZXIgdG8gKiphc3Nlc3MqKiAoYW5kIG1heWJlLCBfcmVtb3ZlXykgdGhpcyBiaWFzLCB3ZSBoYXZlIHRvICoqcXVhbnRpZnkqKiBpdC4KCi0gICBXZSB3aWxsIHBlcmZvcm0gdGhpcyBlc3RpbWF0aW9uIHRoYW5rcyB0byBoZXVyaXN0aWNzIGJhc2VkIG9uICoqa25vd2xlZGdlKiogOiBTZXVyYXQgaW5jbHVkZXMgYSBtZXRob2QgdGhhdCBldmFsdWF0ZXMgdGhlIGNlbGwgY3ljbGUgcGhhc2Ugb2YgY2VsbHMgdGhyb3VnaCAqKnNjb3JlcyBmb3IgdGhlIFMgYW5kIEcyTSBwaGFzZXMqKiwgZWFjaCBiYXNlZCBvbiBwaGFzZS1zcGVjaWZpYyAqKmdlbmUgc2lnbmF0dXJlcyoqLgoKLSAgIEZvciB0aGlzIHN0ZXAsIHdlIHdpbGwgdXNlIGFkZGl0aW9uYWwgZ2VuZSBsaXN0cyBmcm9tIGtub3dsZWRnZQogICAgKGNlbGwgY3ljbGUgcGhhc2UpLCBob3N0ZWQgaW4gYSBaZW5vZG8gcmVzcG9zaXRvcnkgKElkIDoKICAgIFsxNDAzNzM1NV0oaHR0cHM6Ly96ZW5vZG8ub3JnL3JlY29yZHMvMTQwMzczNTUgIlplbm9kbyBnZW5lIGxpc3RzIil7dGFyZ2V0PSJfYmxhbmsifSkKCiMjIERvd25sb2FkIGdlbmUgbGlzdHMKCi0gICBXZSB3aWxsIGRpcmVjdGx5IHJldHJpZXZlIGRhdGEgZnJvbSBaZW5vZG8gdG8geW91ciBgaW5wdXRfZGlyYCA6CgogICAgYGBge3IgZGx6ZW5nbH0KICAgICMgZGx6ZW5nbAogICAgCiAgICAjIyBaZW5vZG8gSUQKICAgIHplbl9pZCA8LSAnMTQxMDE1MDYnCiAgICAKICAgICMjIyBOYW1lZCBmaWxlcyAod2lsbCBiZSB1c2VkIGxhdGVyIG9uICEpCiAgICBjY19maWxlIDwtICJtdXNfbXVzY3VsdXNfU2V1cmF0X2NjLmdlbmVzXzIwMTkxMDMxLnJkcyIKICAgIAogICAgIyMgRmlsZW5hbWUocykgdG8gcmV0cmlldmUKICAgIHRvZ2V0X2ZpbGVzIDwtIGMoY2NfZmlsZSkKICAgIAogICAgIyMgRm9sZGVyIHRvIHN0b3JlIHJldHJpZXZlZCBmaWxlcwogICAgbG9jYWxfZm9sZGVyIDwtIGdsaXN0X2RpcgogICAgCiAgICAjIyBVc2UgbG9jYWwgYmFja3VwID8KICAgIGJhY2t1cCA8LSBGQUxTRQogICAgaWYoYmFja3VwKSBtZXNzYWdlKCJVc2luZyBsb2NhbCBiYWNrdXAgISIpCiAgICAKICAgICMjIEZvcmNlIGRvd25sb2FkID8KICAgIGZvcmNlIDwtIEZBTFNFCiAgICBpZihmb3JjZSkgbWVzc2FnZSgiRm9yY2luZyAocmUpZG93bmxvYWQgISIpCiAgICAKICAgICMjIyBEZWZpbmUgcmVtb3RlIGZvbGRlcgogICAgcmVtb3RlX2ZvbGRlciA8LSBpZiAoYmFja3VwKSB7CiAgICAgICIvc2hhcmVkL3Byb2plY3RzLzI0MjJfZWJhaWlfbjEvYXRlbGllcl9zY3JuYXNlcS9URC9SRVNPVVJDRVMvR0VORUxJU1RTLyIKICAgIH0gIGVsc2UgewogICAgICBwYXN0ZTAoImh0dHBzOi8vemVub2RvLm9yZy9yZWNvcmRzLyIsIHplbl9pZCwgIi9maWxlcy8iKQogICAgfQogICAgCiAgICAjIyMgUmVjb25zdHJ1Y3QgdGhlIGlucHV0IHBhdGhzCiAgICByZW1vdGVfcGF0aCA8LSBwYXN0ZTAocmVtb3RlX2ZvbGRlciwgIi8iLCB0b2dldF9maWxlcykKICAgIAogICAgIyMjIFJlY29uc3RydWN0IHRoZSBvdXRwdXQgcGF0aHMKICAgIGxvY2FsX3BhdGggPC0gcGFzdGUwKGxvY2FsX2ZvbGRlciwgIi8iLCB0b2dldF9maWxlcykKICAgIAogICAgIyMgUmV0cmlldmUgZmlsZXMgKGlmIHRoZXkgZG9uJ3QgZXhpc3QpLCBpbiBsb29wCiAgICBmb3IgKHRnIGluIHNlcV9hbG9uZyh0b2dldF9maWxlcykpIHsKICAgICAgIyMgSWYgdGhlIGZpbGUgZG9lcyBub3QgbG9jYWxseSBleGlzdAogICAgICBpZiAoIWZpbGUuZXhpc3RzKGxvY2FsX3BhdGhbdGddKSB8IGZvcmNlKSB7CiAgICAgICAgIyMgUmV0cmlldmUgZGF0YQogICAgICAgIGlmKGJhY2t1cCkgewogICAgICAgICAgZmlsZS5jb3B5KGZyb20gPSByZW1vdGVfcGF0aFt0Z10sCiAgICAgICAgICAgICAgICAgICAgdG8gPSBsb2NhbF9wYXRoW3RnXSkKICAgICAgICB9IGVsc2UgewogICAgICAgICAgZG93bmxvYWQuZmlsZSh1cmwgPSByZW1vdGVfcGF0aFt0Z10sIAogICAgICAgICAgICAgICAgICAgICAgICBkZXN0ZmlsZSA9IGxvY2FsX3BhdGhbdGddKQogICAgICAgIH0KICAgICAgICAjIyBDaGVjayBpZiBkb3dubG9hZGVkIGZpbGVzIGV4aXN0IGxvY2FsbHkKICAgICAgICBpZihmaWxlLmV4aXN0cyhsb2NhbF9wYXRoW3RnXSkpIG1lc3NhZ2UoIlx0T0siKQogICAgICB9IGVsc2UgbWVzc2FnZShwYXN0ZTAodG9nZXRfZmlsZXNbdGddLCAiIGFscmVhZHkgZG93bmxvYWRlZCAhIikpCiAgICB9CiAgICAKICAgIGBgYAoKPGJyPgoKIyMgTG9hZCBnZW5lIGxpc3RzCgpgYGB7ciBjY19sb2FkfQojIGNjX2xvYWQKCiMjIFRoZSBjZWxsIGN5Y2xlIGdlbmUgbGlzdHMgZmlsZQpjY19maWxlIDwtIHBhc3RlMChnbGlzdF9kaXIsIAogICAgICAgICAgICAgICAgICAiLyIsIAogICAgICAgICAgICAgICAgICBjY19maWxlKQoKIyMgTG9hZCB0aGUgY2VsbCBjeWNsZSByZWZlcmVuY2UgZ2VuZXMgbGlzdHMKY2NfZ2VuZXMgPC0gcmVhZFJEUyhmaWxlID0gY2NfZmlsZSkKCiMjIEhhdmUgYSBsb29rIG9uIGl0cyBjb250ZW50CnN0cihjY19nZW5lcykKCmBgYAoKPGJyPgoKLSAgICoqUXVlc3Rpb24qKgogICAgCiAgICBgYGB7ciBxX2NjX2NoZWNrLCBjbGFzcy5zb3VyY2U9InF1ZXN0aW9uIiwgZXZhbCA9IEZBTFNFfQogICAgIyBxX2NjX2NoZWNrCiAgICAKICAgIEFzIGV4cGxhaW5lZCwgd2Ugd2lsbCB1c2UgZ2VuZSBsaXN0cyBleHRyYWN0ZWQgZnJvbSBjb21tdW5pdHkga25vd2xlZGdlLiAKICAgIE91ciBkYXRhIGNvbnRhaW4gdmFsdWVzIGZvciBnZW5lcyBhbHNvLCBzbyB3ZSB3aWxsIGNyb3NzIHRoZW0uIAogICAgCiAgICBJcyB0aGVyZSBzb21ldGhpbmcgd2Ugc2hvdWxkIGNoZWNrID8KICAgIAogICAgYGBgCiAgICAKICAgIDxicj4KICAgIAogICAgYGBge3IgYV9jY19jaGVjaywgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImFuc3dlciIpLCBldmFsID0gRkFMU0V9CiAgICAjIGFfY2NfY2hlY2sKICAgIAogICAgIyMgLiBXZSBtYXkgY2hlY2sgaWYgdGhlIGdlbmVzIGluIG91ciBnZW5lIGxpc3RzIGFyZSBlZmZlY3RpdmVseSAKICAgICMjICAgcHJlc2VudCBpbiBvdXIgU2V1cmF0IG9iamVjdCAhCiAgICAjIwogICAgIyMgLiBUaGlzIGlzIGV4cGVjdGVkLCBhcyB3ZSByZW1vdmVkIGZlYXR1cmVzIHdpdGggYWxtb3N0IG5vIGV4cHJlc3Npb24KICAgIAogICAgYGBgCiAgICAKICAgIDxicj4KCioqQmV5b25kKiogOgoKMS4gIFdyaXRlIGEgY29kZSB0aGF0IHBlcmZvcm1zIHRoaXMgY2hlY2sKCjIuICBBZGQgYSBjb2RlIHRvIGFkanVzdCB0aGUgY29udGVudCBvZiB0aGUgZ2VuZSBsaXN0cyBhY2NvcmRpbmdseSAocmVtb3ZlIGdlbmVzIGZyb20gdGhlIGdlbmUgbGlzdHMgdGhhdCBhcmUgbm90IHByZXNlbnQgaW4gb3VyIGRhdGFzZXQpCgozLiAgKihOT1RFIDogdGhpcyBpcyBhY3R1YWxseSBub3QgbmVlZGVkIGFzIGFscmVhZHkgY2hlY2tlZCBhbmQKICAgIGNvcnJlY3RlZCBieSB0aGUgY2VsbCBjeWxlIGVzdGltYXRpb24gbWV0aG9kIHdlIHdpbGwgdXNlKSoKCjxicj4KCmBgYHtyIGJfY2NfYW5zd2VyLCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYmV5b25kIiksIGNsYXNzLm91dHB1dD0iYmV5b25kbyIsIGV2YWwgPSBGQUxTRX0KIyBiX2NjX2Fuc3dlcgoKIyMgQ2hlY2sgaWYgb3VyIGRhdGEgZ2VuZXMgY292ZXIgdGhlIGdlbmUgbGlzdHMKbGFwcGx5KGNjX2dlbmVzLCBmdW5jdGlvbih4KSB4ICVpbiUgcm93bmFtZXMoc29iaikpCgojIyBSZW1vdmUgZ2VuZXMgbm90IGluIHNvYmoKY2NfZ2VuZXMgPC0gbGFwcGx5KGNjX2dlbmVzLCBmdW5jdGlvbihnbCkgeyBnbFtnbCAlaW4lIHJvd25hbWVzKHNvYmopXSB9KQoKIyMgQ2hlY2sgdGhlIG1vZGlmaWNhdGlvbgpzdHIoY2NfZ2VuZXMpCgojIyBDaGVjayB0aGF0IGFsbCBnZW5lcyBhcmUgYXZhaWxhYmxlIGluIHNvYmoKYWxsKHVuaXF1ZSh1bmxpc3QoY2NfZ2VuZXMpKSAlaW4lIHJvd25hbWVzKHNvYmopKQoKYGBgCgo8YnI+PGJyPgoKIyMgRXN0aW1hdGlvbgoKLSAgIExldCdzIHBlcmZvcm0gdGhpcyBlc3RpbWF0aW9uLiBCdXQgaG93ID8KCiAgICBgYGB7ciBoX0NlbGxDeWNsZVNjb3JpbmcsIGNsYXNzLnNvdXJjZSA9ICJub3RydW4iLCBldmFsID0gRkFMU0V9CiAgICAjIGhfQ2VsbEN5Y2xlU2NvcmluZwogICAgCiAgICAjIyBSZWFkaW5nIHRoZSBmdW5jdGlvbiBoZWxwIHBhZ2UKICAgID9TZXVyYXQ6OkNlbGxDeWNsZVNjb3JpbmcKICAgIAogICAgYGBgCiAgICAKPGJyPjxicj48YnI+CgotICAgUnVuIHRoZSBjZWxsLWN5Y2xlIGVzdGltYXRpb24KCiAgICBgYGB7ciBjY19ydW59CiAgICAjIGNjX3J1bgogICAgCiAgICAjIyBDcmVhdGluZyBhIHRlbXBvcmFyeSBvYmplY3QKICAgIHNvYmpfQ0NTIDwtIHNvYmoKICAgIAogICAgc29ial9DQ1MgPC0gU2V1cmF0OjpOb3JtYWxpemVEYXRhKG9iamVjdCA9IHNvYmpfQ0NTLCB2ZXJib3NlID0gRkFMU0UpCiAgICAKICAgICMjIFBlcmZvcm0gdGhlIGVzdGltYXRpb24KICAgIHNldC5zZWVkKG15X3NlZWQpCiAgICBzb2JqX0NDUyA8LSBTZXVyYXQ6OkNlbGxDeWNsZVNjb3JpbmcoCiAgICAgIG9iamVjdCA9IHNvYmpfQ0NTLCAKICAgICAgcy5mZWF0dXJlcyA9IGNjX2dlbmVzJHMuZ2VuZXMsIAogICAgICBnMm0uZmVhdHVyZXMgPSBjY19nZW5lcyRnMm0uZ2VuZXMsIAogICAgICAjIGFzc2F5ID0gJ1JOQSwgIyMgU2V1cmF0IHY0CiAgICAgIGxheWVyID0gJ1JOQScsICMjIFNldXJhdCB2NQogICAgICBuYmluID0gMjQsIAogICAgICBzZWVkID0gbXlfc2VlZCkKICAgIAogICAgIyMgVHJhbnNmZXJpbmcgc2NvcmUgZnJvbSB0aGUgdGVtcCBvYmplY3QgdG8gdGhlIHJlYWwgb25lCiAgICBzb2JqJENDX1NldXJhdF9TLlNjb3JlIDwtIHNvYmpfQ0NTJFMuU2NvcmUKICAgIHNvYmokQ0NfU2V1cmF0X0cyTS5TY29yZSA8LSBzb2JqX0NDUyRHMk0uU2NvcmUKICAgIHNvYmokQ0NfU2V1cmF0X1BoYXNlIDwtIGFzLmZhY3Rvcihzb2JqX0NDUyRQaGFzZSkKICAgIAogICAgIyMgQWRkICJTIG1pbnVzIEcyTSIgc2NvcmUKICAgIHNvYmokQ0NfU2V1cmF0X1NtRzJNLlNjb3JlIDwtIHNvYmokQ0NfU2V1cmF0X1MuU2NvcmUgLSBzb2JqJENDX1NldXJhdF9HMk0uU2NvcmUKICAgIAogICAgIyMgRGlzY2FyZCB0ZW1wIG9iamVjdAogICAgcm0oc29ial9DQ1MpCiAgICBgYGAKCjxicj4KCi0gICBEZXNjcmlwdGlvbiBvZiB0aGUgb2JqZWN0IHRvIHNlZSB0aGUgZGF0YSBhZGRlZAoKICAgIGBgYHtyIGRlc2NfY2MsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIiwgZXZhbCA9IEZBTFNFfQogICAgIyBkZXNjX2NjCiAgICAKICAgIFZpZXcoKQogICAgCiAgICBgYGAKCjxicj4KCiMjIFZpc3VhbGl6YXRpb24gCgpBcyB1c3VhbCwgd2UgY2FuIHZpc3VhbGl6ZSB0aGUgcmVzdWx0cyBhcyB2aW9saW5zIDoKCmBgYHtyIHZsbmNjLCBmaWcud2lkdGggPSAxMiwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojdmxuY2MKClNldXJhdDo6VmxuUGxvdChvYmplY3QgPSBzb2JqLAogICAgICAgICAgICAgICAgZmVhdHVyZXMgPSBjKCJDQ19TZXVyYXRfUy5TY29yZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkNDX1NldXJhdF9HMk0uU2NvcmUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICJDQ19TZXVyYXRfU21HMk0uU2NvcmUiKSkKCmBgYAoKPGJyPgoKQnV0IGl0J3Mgbm90IHRoYXQgZWFzeSB0byBpbnRlcnByZXQuLi4gCgpMZXQncyBwbG90IGl0IGluIHRoZSBjZWxsIHNwYWNlCgojIyB7LnVubnVtYmVyZWQgLnRhYnNldCAudGFic2V0LWZhZGUgLnRhYnNldC1waWxsc30KCiMjIyBFc3RpbWF0ZWQgY2VsbCBwaGFzZSA6CgpgYGB7ciB2aXpjY3B9CiMgdml6Y2NwCgojIyBDcmVhdGluZyBhIHRlbXBvcmFyeSBvYmplY3QuCiMjIFRoaXMgdGltZSB3ZSB3aWxsIGtlZXAgaXQgdG8gc3BlZWQgdXAgZnVydGhlciBwbG90cyAobmV4dCBjaHVua3MpCnNvYmpfVklaX2NjIDwtIHNvYmoKCiMjIFVzdWFsIHNldCBvZiAiVSBjYW4ndCBDIG1lIiBmdW5jdGlvbnMKc29ial9WSVpfY2MgPC0gU2V1cmF0OjpOb3JtYWxpemVEYXRhKG9iamVjdCA9IHNvYmpfVklaX2NjLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX2NjIDwtIFNldXJhdDo6U2NhbGVEYXRhKG9iamVjdCA9IHNvYmpfVklaX2NjLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX2NjIDwtIFNldXJhdDo6RmluZFZhcmlhYmxlRmVhdHVyZXMob2JqZWN0ID0gc29ial9WSVpfY2MsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfY2MgPC0gU2V1cmF0OjpSdW5QQ0Eob2JqZWN0ID0gc29ial9WSVpfY2MsIG5wY3MgPSAyMSwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9jYyA8LSBTZXVyYXQ6OlJ1blVNQVAob2JqZWN0ID0gc29ial9WSVpfY2MsIGRpbXMgPSBjKDE6MjApLCB2ZXJib3NlID0gRkFMU0UpCgojIyBDZWxsIGN5Y2xlIHBoYXNlIGZlYXR1cmUgcGxvdApkcF9jY3AgPC0gU2V1cmF0OjpEaW1QbG90KAogIG9iamVjdCA9IHNvYmpfVklaX2NjLCAKICByZWR1Y3Rpb24gPSAndW1hcCcsIAogIGdyb3VwLmJ5ID0gJ0NDX1NldXJhdF9QaGFzZScpICsgU2V1cmF0OjpEYXJrVGhlbWUoKQpwcmludChkcF9jY3ApCgpgYGAKCjxicj4KCiMjIyBTbUcyTSAoIlMgbWludXMgRzJNIikgc2NvcmUgOgoKYGBge3Igdml6Y2NzfQojIHZpemNjcwoKIyMgU21HMk0gZmVhdHVyZSBwbG90CmZwX2NjIDwtIFNldXJhdDo6RmVhdHVyZVBsb3QoCiAgb2JqZWN0ID0gc29ial9WSVpfY2MsIAogIHJlZHVjdGlvbiA9ICd1bWFwJywKICBmZWF0dXJlcyA9ICdDQ19TZXVyYXRfU21HMk0uU2NvcmUnKSArIFNldXJhdDo6RGFya1RoZW1lKCkKcHJpbnQoZnBfY2MpCgojIyBXZSBjYW4gZGlzY2FyZCB0aGUgdGVtcCBvYmplY3QKcm0oc29ial9WSVpfY2MpCgpgYGAKCjxicj4KIyMgey51bm51bWJlcmVkfQoKKipRdWVzdGlvbnMqKgoKYGBge3IgcV9jZWxsMSwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KIyBxX2NlbGwxCgpEb2VzIHRoZSBzdHJ1Y3R1cmUgb2YgdGhlIGNlbGxzIGluIHRoaXMgcmVwcmVzZW50YXRpb24gc2VlbSB0byBoYXZlIAphIGxpbmsgd2l0aCB0aGVzZSBjZWxsIGN5Y2xlIHBoYXNlcy9zY29yZXMgPwoKYGBgCgo8YnI+CgpgYGB7ciBhX2NlbGwxLCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYW5zd2VyIiksIGV2YWwgPSBGQUxTRX0KIyBhX2NlbGwxCgojIyBJdCdzIGFuIGFic29sdXRlIHllcyAhCgpgYGAKCjxicj4KCmBgYHtyIHFfY2VsbDIsIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiMgcV9jZWxsMgoKRG8geW91IHRoaW5rIHRoaXMgaXMgdGhlIHJlc3VsdCBvZiBhbiBhcnRpZmFjdCwgb3IKYmlvbG9neS1yZWxhdGVkID8KCmBgYAoKPGJyPgoKYGBge3IgYV9jZWxsMiwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImFuc3dlciIpLCBldmFsID0gRkFMU0V9CiMgYV9jZWxsMgoKIyMKIyMgICAgwq9cXyjjg4QpXy/CrwojIwoKYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBTYXZlIHRoZSBTZXVyYXQgb2JqZWN0CgpXZSB3aWxsIHNhdmUgb3VyIFNldXJhdCBvYmplY3QgdGhhdCBub3cgY29udGFpbnMgdGhlIGNlbGwgY3ljbGUKc3RhdGVzL3Njb3JlcyA6CgpgYGB7ciBzYXZlcmRzMiwgZm9sZC5vdXRwdXQgPSBGQUxTRX0KIyBzYXZlcmRzMgoKIyMgU2F2ZSBvdXIgU2V1cmF0IG9iamVjdCAocmljaCBuYW1pbmcpCm91dF9uYW1lIDwtIHBhc3RlMCgKICAgICAgICAgIG91dHB1dF9kaXIsICIvIiwgcGFzdGUoCiAgICAgICAgICAgIGMoIjA0IiwgU2V1cmF0OjpQcm9qZWN0KHNvYmopLCAiUzUiLCAKICAgICAgICAgICAgICAiQ0MiLCBwYXN0ZSgKICAgICAgICAgICAgICAgIGRpbShzb2JqKSwgCiAgICAgICAgICAgICAgICBjb2xsYXBzZSA9ICcuJwogICAgICAgICAgICAgICkKICAgICAgICAgICAgKSwgY29sbGFwc2UgPSAiXyIpLAogICAgICAgICAgICAiLlJEUyIpCgojIyBDaGVjawpwcmludChvdXRfbmFtZSkKCiMjIFdyaXRlIG9uIGRpc2sKc2F2ZVJEUyhvYmplY3QgPSBzb2JqLCAKICAgICAgICBmaWxlID0gb3V0X25hbWUpCgpgYGAKCjxicj4KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgo8YnI+CgojIENlbGwgZG91YmxldHMKCldlIHdpbGwgdXNlICoqdHdvIGRpZmZlcmVudCBtZXRob2RzKiogdG8gZGV0ZWN0IGFuZCByZW1vdmUgY2VsbCBkb3VibGV0cyA6CgotICAgW2BzY2RzYF0oaHR0cHM6Ly93d3cuYmlvY29uZHVjdG9yLm9yZy9wYWNrYWdlcy9yZWxlYXNlL2Jpb2MvaHRtbC9zY2RzLmh0bWwpe3RhcmdldD0iX2JsYW5rIn0gKGluIGl0cyAiaHlicmlkIiBtb2RlKSA6IG1vcmUgZWZmaWNpZW50IGF0IGRldGVjdGluZyAqKmhvbW90eXBpYyoqIGRvdWJsZXRzCgotICAgW2BzY0RibEZpbmRlcmBdKGh0dHBzOi8vd3d3LmJpb2NvbmR1Y3Rvci5vcmcvcGFja2FnZXMvcmVsZWFzZS9iaW9jL2h0bWwvc2NEYmxGaW5kZXIuaHRtbCl7dGFyZ2V0PSJfYmxhbmsifSA6IG1vcmUgZWZmaWNpZW50IGF0IGRldGVjdGluZyAqKmhldGVyb3R5cGljKiogZG91YmxldHMKCk5vbmUgb2YgdGhlIG1ldGhvZHMgYWNjZXB0cyBhIGBTZXVyYXRPYmplY3RgIGFzIGlucHV0LCBidXQgYSBbYFNpbmdsZUNlbGxFeHBlcmltZW50YF0oaHR0cHM6Ly93d3cuYmlvY29uZHVjdG9yLm9yZy9wYWNrYWdlcy9yZWxlYXNlL2Jpb2MvaHRtbC9TaW5nbGVDZWxsRXhwZXJpbWVudC5odG1sKXt0YXJnZXQ9Il9ibGFuayJ9IG9iamVjdC4KCkhvcGVmdWxseSA6CgotICAgU2V1cmF0IGhhcyBhIGZ1bmN0aW9uIHRvIHBlcmZvcm0gdGhlIGNvbnZlcnNpb24KLSAgIFRoZSBvdXRwdXQgcmVzdWx0cyBjYW4gYmUgaW50ZWdyYXRlZCBpbnRvIG91ciBTZXVyYXQgb2JqZWN0IHdpdGggZWFzZQoKIyMgVHdvIG1ldGhvZHMgey50YWJzZXQgLnRhYnNldC1mYWRlIC50YWJzZXQtcGlsbHN9CgojIyMgYHNjZHNgCgpgYGB7ciBzY2RzfQojIHNjZHMKCiMjIEZpeCBzZWVkCnNldC5zZWVkKG15X3NlZWQpCgojIyBSdW4gc2Nkcwpzb2JqJGRvdWJsZXRfc2Nkcy5oeWJyaWQgPC0gdW5uYW1lKAogIHNjZHM6OmN4ZHNfYmNkc19oeWJyaWQoCiAgICBTZXVyYXQ6OmFzLlNpbmdsZUNlbGxFeHBlcmltZW50KAogICAgICBzb2JqLCBhc3NheSA9ICJSTkEiKSkkaHlicmlkX3Njb3JlID4gMSkKCiMjIENvbnRpbmdlbmNpZXMKdGFibGUoc29iaiRkb3VibGV0X3NjZHMuaHlicmlkKQoKYGBgCgoKIyMjIGBzY0RibEZpbmRlcmAKCmBgYHtyIHNjZGJsfQojIHNjZGJsCgojIyBGaXggc2VlZApzZXQuc2VlZChteV9zZWVkKQoKIyMgUnVuIHNjRGJsRmluZGVyICh3aGljaCBuZWVkcyBhbm90aGVyIG9iamVjdCB0eXBlKQpzb2JqJGRvdWJsZXRfc2NEYmxGaW5kZXIgPC0gc2NEYmxGaW5kZXI6OnNjRGJsRmluZGVyKAogIHNjZSA9IFNldXJhdDo6YXMuU2luZ2xlQ2VsbEV4cGVyaW1lbnQoCiAgICB4ID0gc29iaiwgCiAgICBhc3NheSA9ICJSTkEiKSwKICByZXR1cm5UeXBlID0gInRhYmxlIikkY2xhc3MgPT0gImRvdWJsZXQiCgojIyBDb250aW5nZW5jaWVzCnRhYmxlKHNvYmokZG91YmxldF9zY0RibEZpbmRlcikKCmBgYAoKIyAgey51bm51bWJlcmVkfQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCjxicj4KCiMjIE1lcmdlIHJlc3VsdHMKCldlIG1lcmdlIHJlc3VsdHMgb2YgdGhlIHR3byBtZXRob2RzCgpgYGB7ciBkYmxtZXJnZX0KIyBkYmxtZXJnZQoKIyMgTG9naWNhbCB1bmlvbiBvZiBib3RoIG1ldGhvZHMKc29iaiRkb3VibGV0X3VuaW9uIDwtIHNvYmokZG91YmxldF9zY2RzLmh5YnJpZCB8IHNvYmokZG91YmxldF9zY0RibEZpbmRlcgoKIyMgUXVhbnRpZnkgZG91YmxldHMKdGFibGUoc29iaiRkb3VibGV0X3VuaW9uKQoKYGBgCgo8YnI+CgpXZSBjYW4gYXNzZXNzICoqdG9vbC1zcGVjaWZpYyoqIGFuZCAqKmNvbW1vbioqIGRvdWJsZXRzCgpgYGB7ciBkYmxfdHlwZXMsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIn0KIyBkYmxfdHlwZXMKCiMjIyBTaW5nbGV0cyBieSBkZWZhdWx0CnNvYmokZG91YmxldF92aXogPC0gInNpbmdsZXQiCgojIyMgVW5pb24Kc29iaiRkb3VibGV0X3Zpeltzb2JqJGRvdWJsZXRfdW5pb25dIDwtICJib3RoIgoKIyMjIHNjZHMtc3BlY2lmaWMKc29iaiRkb3VibGV0X3Zpeltzb2JqJGRvdWJsZXRfc2Nkcy5oeWJyaWQgJiAhc29iaiRkb3VibGV0X3NjRGJsRmluZGVyXSA8LSAic2NkcyIKCiMjIyBzY0RibEZpbmRlci1zcGVjaWZpYwpzb2JqJGRvdWJsZXRfdml6W3NvYmokZG91YmxldF9zY0RibEZpbmRlciAmICFzb2JqJGRvdWJsZXRfc2Nkcy5oeWJyaWRdIDwtICJzY0RibEZpbmRlciIKCiMjIENvbnZlcnQgdG8gZmFjdG9yCnNvYmokZG91YmxldF92aXogPC0gYXMuZmFjdG9yKHNvYmokZG91YmxldF92aXopCgojIyBDb250aW5nZW5jaWVzCnRhYmxlKHNvYmokZG91YmxldF92aXopCgpgYGAKCjxicj4KCioqQmV5b25kKiogOiBDcmVhdGUgYW4gdXBzZXQtcGxvdCBmb3IgdGhlIGRvdWJsZXQgc3RhdHVzIGFjY29yZGluZyB0byB0aGUgdHdvIG1ldGhvZHMgdXNlZAoKYGBge3IgYl91cHNldCwgY2xhc3Muc291cmNlPWMoImZvbGQtaGlkZSIsICJiZXlvbmQiKSwgY2xhc3Mub3V0cHV0PSJiZXlvbmRvIn0KIyBiX3Vwc2V0CgojIyBCdWlsZCBhIGxpc3Qgb2YgY2VsbHMgdGFnZ2VkIGJ5IG9uZSB0b29sLCB0aGUgb3RoZXIsIG9yIGJvdGgKZGJsX2xpc3QgPC1saXN0KAogICJzY2RzIiA9IGNvbG5hbWVzKHNvYmopW3NvYmokZG91YmxldF9zY2RzLmh5YnJpZCAmICFzb2JqJGRvdWJsZXRfc2NEYmxGaW5kZXJdLAogICJzY0RibEZpbmRlciIgPSBjb2xuYW1lcyhzb2JqKVshc29iaiRkb3VibGV0X3NjZHMuaHlicmlkICYgc29iaiRkb3VibGV0X3NjRGJsRmluZGVyXSwKICAiYm90aCIgPSBjb2xuYW1lcyhzb2JqKVtzb2JqJGRvdWJsZXRfc2Nkcy5oeWJyaWQgJiBzb2JqJGRvdWJsZXRfc2NEYmxGaW5kZXJdKQoKIyMgRHJhdyB0aGUgdXBzZXQgcGxvdApVcFNldFI6OnVwc2V0KGRhdGEgPSBVcFNldFI6OmZyb21MaXN0KGRibF9saXN0KSwgCiAgICAgICAgICAgICAgbmludGVyc2VjdHMgPSBOQSwgCiAgICAgICAgICAgICAgc2V0cyA9IHJldihuYW1lcyhkYmxfbGlzdCkpLAogICAgICAgICAgICAgIGtlZXAub3JkZXIgPSBUUlVFLAogICAgICAgICAgICAgIG9yZGVyLmJ5ID0gImZyZXEiKQoKYGBgCgo8YnI+CgpOb3cgd2UgY2FuIHJlbW92ZSBiYXJjb2RlcyBpZGVudGlmaWVkIGFzIGNlbGwgZG91YmxldHMsIGFuZCB2aXN1YWxpemUgdGhlIGNlbGwgc3BhY2UgYmVmb3JlIGFuZCBhZnRlci4KCiMjIERvdWJsZXRzIGZpbHRlcmluZyB7LnRhYnNldCAudGFic2V0LWZhZGUgLnRhYnNldC1waWxsc30KCiMjIyBCRUZPUkUKCmBgYHtyIGRibHZpejF9CiMgZGJsdml6MQoKIyMjIERvdWJsZXRzIHZpeiAoYmVmb3JlIHJlbW92YWwpCgojIyBDcmVhdGUgYSB0ZW1wb3Jhcnkgb2JqZWN0CnNvYmpfVklaX2RibCA8LSBzb2JqCgojIyBVc3VhbCBzZXQgb2YgIlUgY2FuJ3QgQyBtZSIgZnVuY3Rpb25zCnNvYmpfVklaX2RibCA8LSBTZXVyYXQ6Ok5vcm1hbGl6ZURhdGEob2JqZWN0ID0gc29ial9WSVpfZGJsLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX2RibCA8LSBTZXVyYXQ6OlNjYWxlRGF0YShvYmplY3QgPSBzb2JqX1ZJWl9kYmwsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfZGJsIDwtIFNldXJhdDo6RmluZFZhcmlhYmxlRmVhdHVyZXMob2JqZWN0ID0gc29ial9WSVpfZGJsLCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfVklaX2RibCA8LSBTZXVyYXQ6OlJ1blBDQShvYmplY3QgPSBzb2JqX1ZJWl9kYmwsIG5wY3MgPSAyMSwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9kYmwgPC0gU2V1cmF0OjpSdW5VTUFQKG9iamVjdCA9IHNvYmpfVklaX2RibCwgZGltcyA9IGMoMToyMCksIHZlcmJvc2UgPSBGQUxTRSkKCiMjIERvdWJsZXRzIHZpc3VhbGl6YXRpb24KZHBfZGJsIDwtIFNldXJhdDo6RGltUGxvdCgKICBvYmplY3QgPSBzb2JqX1ZJWl9kYmwsIAogIHJlZHVjdGlvbiA9ICd1bWFwJywKICBncm91cC5ieSA9ICdkb3VibGV0X3ZpeicpICsgU2V1cmF0OjpEYXJrVGhlbWUoKQoKIyMgUGxvdApwcmludChkcF9kYmwpCgojIyBEaXNjYXJkIHRlbXAgb2JqZWN0CnJtKHNvYmpfVklaX2RibCkKCmBgYAoKIyMjIFJlbW92ZSBkb3VibGV0cwoKYGBge3IgZGJscm0sIGZvbGQub3V0cHV0ID0gRkFMU0V9CiMgZGJscm0KCiMjIERpbWVuc2lvbnMgYmVmb3JlIHJlbW92YWwKZGltKHNvYmopCgojIyBQZXJmb3JtIHRoZSByZW1vdmFsCiMgc29iaiA8LSBzb2JqWywhc29iaiRkb3VibGV0X3VuaW9uXQpzb2JqIDwtIHN1YnNldCgKICB4ID0gc29iaiwgCiAgY2VsbHMgPSBjb2xuYW1lcyhzb2JqKVtzb2JqJGRvdWJsZXRfdW5pb25dLCAKICBpbnZlcnQgPSBUUlVFKQoKIyMgRGltZW5zaW9ucyBhZnRlcgpkaW0oc29iaikKCmBgYAoKIyMjIEFGVEVSCgpgYGB7ciBkYmx2aXoyLCBjbGFzcy5zb3VyY2UgPSAibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIn0KIyBkYmx2aXoyCgojIyMgRG91YmxldHMgdml6IChhZnRlciByZW1vdmFsKQoKIyMgQ3JlYXRlIGEgdGVtcG9yYXJ5IG9iamVjdApzb2JqX1ZJWl9kYmxmIDwtIHNvYmoKCiMjIFVzdWFsIHNldCBvZiAiVSBjYW4ndCBDIG1lIiBmdW5jdGlvbnMKc29ial9WSVpfZGJsZiA8LSBTZXVyYXQ6Ok5vcm1hbGl6ZURhdGEob2JqZWN0ID0gc29ial9WSVpfZGJsZiwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX1ZJWl9kYmxmIDwtIFNldXJhdDo6U2NhbGVEYXRhKG9iamVjdCA9IHNvYmpfVklaX2RibGYsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfZGJsZiA8LSBTZXVyYXQ6OkZpbmRWYXJpYWJsZUZlYXR1cmVzKG9iamVjdCA9IHNvYmpfVklaX2RibGYsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfZGJsZiA8LSBTZXVyYXQ6OlJ1blBDQShvYmplY3QgPSBzb2JqX1ZJWl9kYmxmLCBucGNzID0gMjEsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9WSVpfZGJsZiA8LSBTZXVyYXQ6OlJ1blVNQVAob2JqZWN0ID0gc29ial9WSVpfZGJsZiwgZGltcyA9IGMoMToyMCksIHZlcmJvc2UgPSBGQUxTRSkKCiMjIERvdWJsZXRzLWZpbHRlcmVkIHZpc3VhbGl6YXRpb24KZHBfZGJsZiA8LSBTZXVyYXQ6OkRpbVBsb3QoCiAgb2JqZWN0ID0gc29ial9WSVpfZGJsZiwgCiAgcmVkdWN0aW9uID0gJ3VtYXAnLAogIGdyb3VwLmJ5ID0gJ2RvdWJsZXRfdml6JykgKyBTZXVyYXQ6OkRhcmtUaGVtZSgpCgojIyBQbG90CnByaW50KGRwX2RibGYpCgojIyBEaXNjYXJkIHRlbXAgb2JqZWN0CnJtKHNvYmpfVklaX2RibGYpCgpgYGAKCiMjICB7LnVubnVtYmVyZWR9CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKPGJyPgoKTWVyZ2UgcGxvdHMgZm9yIGVhc2Ugb2YgdXNlIDoKCmBgYHtyIGRibF91bWFwcywgZmlnLndpZHRoID0gMTIsIGZpZy5oZWlnaHQgPSA2LCBjbGFzcy5zb3VyY2UgPSAibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIn0KIyBkYmxfdW1hcHMKCiMjIFVzaW5nIHRoZSBwYXRjaHdvcmsgcGFja2FnZSB0byBtZXJnZSBwbG90cyAoYW5kIGdncGxvdDIgdG8gYWRkIHRpdGxlcykKcGF0Y2h3b3JrOjp3cmFwX3Bsb3RzKAogIGxpc3QoCiAgICBkcF9kYmwgJiBnZ3Bsb3QyOjpnZ3RpdGxlKGxhYmVsID0gIkNlbGwgZG91YmxldHMgKHVuZmlsdGVyZWQpIiksCiAgICBkcF9kYmxmICYgZ2dwbG90Mjo6Z2d0aXRsZShsYWJlbCA9ICJDZWxsIGRvdWJsZXRzIChmaWx0ZXJlZCkiKSksIAogIG5yb3cgPSAxKQoKYGBgCgo8YnI+CgoqKlF1ZXN0aW9uKioKCmBgYHtyIHFfZGJsY29tcCwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWw9RkFMU0V9CiMgcV9kYmxjb21wCgpXaGF0IGRvIHlvdSBvYnNlcnZlIHdoZW4gY29tcGFyaW5nIGJlZm9yZSBhbmQgYWZ0ZXIgdGhlIGRvdWJsZXRzIGZpbHRlcmluZyA/CgpgYGAKCmBgYHtyIGFfZGJsY29tcCwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImFuc3dlciIpLCBldmFsPUZBTFNFfQojIGFfZGJsY29tcAoKIyMgLiBBY3R1YWxseSBub3QgdGhhdCBtdWNoICEKIyMKIyMgLiBUaGlzIGlzIGR1ZSB0byB0aGUgZmFjdCB0aGF0IHRoaXMgZGF0YXNldCBkaWQgbm90IGNvbnRhaW4KIyMgICBhIGxvdCBvZiBoZXRlcm90eXBpYyBkb3VibGV0cwojIwojIyAuIFdoZW4gYSBkYXRhc2V0IGlzIGhlYXZpbHkgY29udGFtaW5hdGVkIHdpdGggbG90cyBvZiBoZXRlcm90eXBpYwojIyAgIGRvdWJsZXRzLCB0aGVpciByZW1vdmFsIGNhbiBkcmFzdGljYWxseSByZW1vZGVsIHRoZSB0b3BvbG9neS4KCmBgYAoKPGJyPgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCgojIFNhdmUgdGhlIFNldXJhdCBvYmplY3QKCldlIHdpbGwgc2F2ZSBvdXIgU2V1cmF0IG9iamVjdCB0aGF0IGlzIG5vdyBmaWx0ZXJlZCBmb3IgZG91YmxldHMgOgoKYGBge3Igc2F2ZXJkczN9CiMgc2F2ZXJkczMKCiMjIFNhdmUgb3VyIFNldXJhdCBvYmplY3QgKHJpY2ggbmFtaW5nKQpvdXRfbmFtZSA8LSBwYXN0ZTAoCiAgICAgICAgICBvdXRwdXRfZGlyLCAiLyIsIHBhc3RlKAogICAgICAgICAgICBjKCIwNSIsIFNldXJhdDo6UHJvamVjdChzb2JqKSwgIlM1IiwgCiAgICAgICAgICAgICAgIkRvdWJsZXRzLkZpbHRlcmVkIiwgcGFzdGUoCiAgICAgICAgICAgICAgICBkaW0oc29iaiksIAogICAgICAgICAgICAgICAgY29sbGFwc2UgPSAnLicKICAgICAgICAgICAgICApCiAgICAgICAgICAgICksIGNvbGxhcHNlID0gIl8iKSwKICAgICAgICAgICAgIi5SRFMiKQoKIyMgQ2hlY2sKcHJpbnQob3V0X25hbWUpCgojIyBXcml0ZSBvbiBkaXNrCnNhdmVSRFMob2JqZWN0ID0gc29iaiwgCiAgICAgICAgZmlsZSA9IG91dF9uYW1lKQoKYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKPGJyPjxicj48YnI+CgojIFJzZXNzaW9uCgpgYGB7ciByc2Vzc2lvbiwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojIHJzZXNzaW9uCgp1dGlsczo6c2Vzc2lvbkluZm8oKQoKYGBgCg==